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ABSTRACT 

Velocity and temperature fluctuation measurements were 
made over the open ocean from instruments mounted on the 
R/V ACANI A . These data were examined to determine the 

validity of present formulations and prediction techniques. 

2 

Values of momentum flux, u* , were inferred from the 
rate of dissipation of turbulent kinetic energy, e. Dissi- 
pation values were obtained from spectral estimates and 
inner scale estimates. Values of u* were examined for 
representativeness on the basis of the cons t ant- f lux assump- 
tion and by comparisons with other studies. The vertical 
variation of the dissipation rate was examined for possible 
effects of stability and wind-wave coupling. 

The momentum flux, computed from spectra, supported the 
cons tant- f lux assumption for neutral conditions. For 
periods of instability, the stability corrections applied 
to the vertical variation of e resulted in the proper ad- 
justments toward the predicted slope. The reductions in 
momentum transfer during periods of stable stratifications 
were consistent with wind-wave coupling effects described 
by Davidson. The shape of a spectrum of temperature fluc- 
tuations was in agreement with predictions. 
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I 



INTRODUCTION 



The structure of turbulence in the marine boundary 
layer has become of increasing importance in naval systems 
design. To gain a better understanding of turbulence sta- 
tistics, such as the spectral distribution of wind fluc- 
tuations, the effects of wind-wave coupling and the 
accuracy of .present prediction expressions, personnel at 
the Naval Postgraduate School, Monterey, California have 
undertaken a study to obtain the shipboard observational 
data necessary for an accurate examination of these prob- 
lems. The group consists of personnel from several depart- 
ments; Meteorology, Mechanical Engineering, Oceanography 
and Physics, whose primary objective is to obtain observa- 
tional data and then to evaluate and improve existing pre- 
diction techniques . 

Several observational experiments have been conducted 
and were followed by data reduction and analyses periods. 
The second observational experiment is the one considered 
in this study and was conducted from 18 September 1973 to 
22 September 1973. Meteorological measurements were made 
from the oceanographic research ship, the R/V ACANIA shown 
in Figure 1, in the vicinity of San Nicolas Island. Fig- 
ure 2 illustrates the "near" ocean environment and the 
ACAN I A 1 s anchorage where the majority of measurements were 
taken. 
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* 




Figure 1. The R/V ACANIA is operated by the Department 
of Oceanography, NPS. 
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Figure 2. San Nicolas Island and ACANIA 1 



Anchorage. 
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The meteorological measurements were of three differ- 
ent types. Two of these pertained to different scales of 
measurements and the third was of the influence of the 
moving platform on mean and turbulent measurements. The 
scales considered were the turbulent scale and that which 
is associated with the mean wind and thermal stratifica- 
tion. 



The scope of this study is confined to turbulent fluc- 
tuations of velocity and temperature, although reference 
will be made to mean wind and thermal stratification mea- 
surements. The latter are important since they are re- 
quired to relate to results from other experiments or for 
extrapolating to regions or conditions where turbulence 
measurements are not made c A nearly unique problem facing 
the Navy is that present expressions for the turbulent 
processes are empirical in nature and have only been pre- 
viously validated, or at least tested, overland. 

The purpose of this investigation is to test several 
expressions in the "near" ocean environment, interpret the 
results and evaluate existing turbulence expressions. To 
accomplish this, it was decided to concentrate the study 
on specific formulations. The formulations include (1) an 
evaluation of the constant flux layer (2) the variations 
of £ (the viscous dissipation of turbulent kinetic energy) 
with height and, where data permitted, (3) the variation 
of temperature fluctuations. 
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II 



OBSERVATIONAL AND THEORETICAL CONS IDERATIONS 



A. GENERAL BACKGROUND 

This study represents one of the first attempts to 
make simultaneous, multiple-level observations over water 
from a ship D Previous measurements have been primarily 
overland or over water from near stable platforms. A 
study which closely approximates this experiment was de- 
scribed by Stegen, et al (1973). The measurements re- 
ported in that study were made from the Scripps Floating 
Instrument Platform (FLIP) during the Barbados Oceano- 
graphic and Meteorological Experiment (BOMEX). From FLIP, 
velocity and temperature measurements were made with a 
roving probe which cycled continuously between measurement 
stations located 2, 4, 7 and 12 meters above mean sea 

level. Results in the present study will be compared to 
those obtained by Stegen, et al, as well as to other con- 
temporary studies . 

The following paragraphs will provide a review of the 
theoretical considerations needed to interpret the data. 
From these considerations, the friction velocity (u*) will 
be estimated from values of the dissipation of turbulent 
kinetic energy (c) and e will be obtained from spectra. 

For this reason we must examine both the turbulent kinetic 
energy balance expression and Kolmogorov's hypotheses. 
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B. TURBULENT BALANCE EXPRESSIONS 

A general expression for the turbulent kinetic energy 
balance has the form 




-V . V . 
i 3 



3u . 
1 

3x . 
3 
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If all eddy transports except those along the vertical are 
neglected, the turbulent kinetic energy balance can be ex- 
pressed as 





9u , 
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( 2 ) 



The time variation of turbulent kinetic energy is thus seen 
to depend on four processes represented by the four terms 
on the right hand side of Equation (2) . The first repre- 
sents production of turbulent kinetic energy by the inter- 
action of the Reynolds stress with the mean shear flow. The 
second represents turbulent kinetic energy released or 
gained due to the work of vertical motions against the 
effects of buoyancy. The third represents viscous dissipa- 
tion of turbulent kinetic energy. The final term represents 
both pressure work and the flux divergence of turbulent 
kinetic energy due to transfer by turbulent motion. 
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To obtain those expressions normally used for estimat- 



ing the shear stress or eddy flux (-V^V^ = ) on the 

basis of turbulent kinetic energy balance, the following 
assumptions and definitions are made: 



i) assume that the turbulent kinetic energy is in 

a statistically steady state 



» v - 

Jt { ~2~ ) = ° 



ii) conditions are horizontally homogeneous for the 
irregularity fields (i.e., spatial variations 
of mean fluctuations are taken to be negligible 
over distances comparable to the largest scales 
of interest in the analysis). 





0 



iii) conditions are near neutral 




o 



iv) define the coordinate system so 



u 



2 



= u . 



0 



The above assumptions and definitions lead to the follow- 
ing turbulent energy balance expression 
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Furthermore, if the divergence terms are negligible, 
Equation (3) becomes 



3 u 



-V, V 



1 3 3x. 



= e 



(4) 
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which states that mechanical production of turbulent 
kinetic energy equals dissipation of turbulent kinetic 
energy . 

Equation (4) can be rewritten as 



u 



2 du. 
9x. 



= e 



(5) 
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where u = -V V . For neutral conditions in a constant- 
* 13 

stress boundary layer the following is an accepted expres- 
sion for the wind shear 



9u _ * 

9x^ kz 



(6) 



where k is Von Karman’s constant and z is the height. 
Combining Equations (5) and (6) yields the following esti- 
mate for u* based on measurements of e at known height z. 



u 



* 



(ezk) 



1/3 



(7) 



Equation (7) can be used to both estimate the value of the 

2 

momentum flux, u* , and to examine if it is constant with 
respect to height, the cons tant-f lux assumption. 

Table I provides the momentum flux calculations pre- 
sented by Stegen, et al with e values derived from time 
derivatives of the velocity fluctuations. These values 
will be compared to those obtained from shipboard measure- 
ments to establish the validity of the expressions and 
assumptions used in this study. 
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TABLE I 



Results of Momentum-Flux Calculations 



1 





Height 

(m) 


kze 

( cni^sec'^ ) 


U * 

1 

(cm sec"- 1 -) 


u* (corrected ) 
(cm sec -1 ) 





3 






11.69 


4.02x10 


15.9 


14.0 


6.70 


6.86xl0 3 


19.0 


17.3 


3.89 


6.36xl0 3 


18.5 


17.4 


2.39 


6.03xl0 3 


18.2 


17.4 



‘'"Results taken from Stegen, et al (1973) 
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Equation (7) can be rewritten as 



Ln G = -Ln z + Ln u* /k 



( 8 ) 



Therefore a plot of Ln z versus Ln e should yield a 
straight line with -1 slope. The intercept of this line 
with the Ln z axis, z = 1, will yield the estimate of 
u^, A plot of this shown by Stegen, et al appears in 
Figure 3. Deviations from this particular relation could 
be interpreted for possible deviation from the constant- 
flux assumption or for influence of wind-wave coupling 
or non-neutral stratifications. 



C. SIMILARITY EXPRESSIONS 

The parameter e can be estimated on the basis of 
similarity expressions for the inertial subrange. Accord- 
ing to Kolmogorov's second hypothesis, the inertial sub- 
range is a region (in wave number space) where there is 
negligible dissipation of energy and where inertial energy 
transfer to higher wave numbers is the dominate feature. 

Figure 4 from Lumley and Panofsky (1964) depicts the 
inertial subrange considering spectral energy transfers. 

Based on a dimensional analysis Kolmogorov's second 
hypothesis yields the following functional relation 



E (k) = =e 2 / V 5/3 



where E(k) represents the three-dimensional energy density, 
k is the wave number and c is an empirical constant. 
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Measured Values of e vs z, the Height Above 
the Mean Sea Surface (after Stegen, et al). 
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Figure 4. A Schematic Drawing Showing Spectral Transfer 
•and Illustrating (A) Energy Containing Region 
(B) Inertial Subrange (C) Dissipation Range 
and (D) Universal Equilibrium Range. 



Unfortunately, E(k) is impossible to measure. It can be 
shown, however, that the one-dimensional energy density, 
$11 (k), which can be measured, also has the form 

<J* 11 (k) = ae 2/3 k" 5/3 (9) 

where a = constant. 

Furthermore, measurements at the scales of interest 
have to be made at a fixed point in the flow and the spec- 
tra therefore are obtained for temporal frequency, f. In 
order to use Equation (9), time and space scales must be 
related using Taylor's (1938) "frozen turbulence" hypo- 
thesis 

k = 2 7T f / U 

where U is the mean wind speed, measured by a cup anemo- 
meter. The term "frozen turbulence' implies that the 
turbulence remains unchanged during the time required for 
it to sweep past the probe. 

The following form of Equation (9) is often utilized 
and will be here, 

f ^ i i ( f ) = k<J) n (k) = ae 2/3 k -2/3 (10) 

The parameter e will be obtained from the spectra using 

Equation (10). The C values obtained from simultaneous, 

multiple level measurements will then be used to estimate 

2 

the value of the momentum flux, u* , from Equation (7). 
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There is another way of estimating £ based on the 
scale associated with the upper limit of the inertial sub- 
range, If the spectrum is in statistical equilibrium it 
is independent of the source of the turbulence and depends 
only on e f y and k. A dimensional analysis yields the 
following expression for the inner scale, 1 q 

1 Q = (Y 3 /£) 1/4 (ID 

where k g ( = 27rl Q /u) is the dissipation wave number and y 
is the kinematic viscosity of air. The dissipation wave 
number can be obtained from spectra which show scale 
separation between pure inertial transfer (-5/3 regions) 
and the viscous dissipation region. Dissipation of turbu- 
lent kinetic energy can then be obtained from Equation (11) 
and compared to those obtained from Equation (10) to estab- 
lish greater confidence in the spectral estimate. 

D. FURTHER THEORETICAL CONSIDERATIONS 

Although the expressions previously described are based 
on accepted and consistent theories, there are several 
serious assumptions which must be considered in applying 
them to observational data. In many instances the theory 
is empirical and the formulating analyses were based on 
overland data. A recent effort in turbulence theory has 
been to prove or disprove the validity of overland formu- 
lations for non-neutral conditions or over water regimes 
where the wave influence may be important. The latter has 
been described by Davidson (1974). 
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In developing the turbulent kinetic energy balance, 
given by Equation (3), the existence of neutral conditions 
was a critical assumption 0 For non-neutral stratifica- 
tions, Equation (7) would include a buoyancy term and be 
of the form 



e 







112 ) 



or 



e - sf- V Ri) 



where L = 
flux , 



- Tu /gkV 0 = Monin-Obuhkov Length; V 0 = heat 



Ri = Richardson Number, 

2 

<f) = empirical function of — , 

1 L 

(f)^ = empirical function of Ri . 

For the non-neutral case, the expression corresponding to 
Equation (8) becomes 



Ln e 



-Ln z + Ln ({^(Ri) + Ln 




(13) 



Possible deviations from the constant flux assumption or 
wind wave coupling influence exhibited by spectral results 
could be examined by testing overland expressions for (J)^(Ri) 
A stability correction could then be applied for periods of 
non-neutral stratifications. A possible effect of the 
instability on the predicted plots of Ln C versus Ln z is 
illustrated in Figure 5. 



2 3 




LnZ 



LnC 




LnZ 



Figure 5. An Illustration of a Possible Effect of 

Stability on the Variation of c with Height. 
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OBSERVATIONS 



A. THE EXPERIMENT 

The data considered in this study were obtained from 
measurements made aboard the R/V AC AN I A , which is operated 
by the Department of Oceanography of the United States 
Naval Postgraduate School. As previously stated, San 
Nicolas Island was the experimental site and the majority 
of data were collected while the ship was at anchorage near 
the northwest tip of the island. Additional data were col- 
lected while the ship was enroute between Monterey and San 
Nicolas Island and also while the ship was transiting to 
the west side of the island. Since the data collected from 
anchorage comprise the largest percentage of successful 
runs and represents all stability conditions, data runs 
selected for analysis in this study are primarily anchorage 
data. Table II is a listing of the various data runs ex- 
amined. Figure 6 shows the ship's track in the vicinity of 
San Nicolas Island and Figure 7 shows the route to and from 
the island. 

B. DESCRIPTION OF EQUIPMENT AND MEASUREMENTS 

Figure 8 shows the sensor placement at various levels 
along with the variables measured at each level. 

1 . Velocity Measurements 

Velocity fluctuation measurements were made at three 
levels on the forward mast and one level on a single bow 
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TABLE II 



Data Periods Considered in This Study 



Date 



Time Digital Channel Parameter 



20 Sept. 73 1103-1138 



20 Sept. 73 2008-2018 



20 Sept. 73 2023-2053 



1 

2 

3 

4 

1 

2 

3 

4 

5 

6 

7 

8 



1 

2 

3 

4 




upper level 
lateral accelera- 
tion 

upper level ver- 
tical acceleration 

upper level ver- 
tical acceleration 
(integrated ) 



u • 



u 



u ' 



The subscript refers to the level at which the 
measurement was made. Level 1 is the uppermost 
level . 
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TABLE II (continued) 



Date 


Time 


Digital Channel 


Parameter 


20 Sept. 73 


2323-2341 


1 


u i 






2 


U 2 






3 


U 3 




* 


4 


upper level fore-aft 
acceleration 






5 


upper level lateral 
acceleration 






6 


T 4 






7 


upper level vertical 
acceleration 








(integrated) 






8 


upper level vertical 
acceleration 


21 Sept. 73 


0637-0652 


1 


upper level fore-aft 
acceleration 






2 


U 3 






3 


u i 






4 


T 4 






5 


t b 






6 


upper level vertical 
acceleration 






7 


upper level lateral 
acceleration 






8 


upper level lateral 
acceleration 








(integrated) 
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Figure 6. ACAHIA's Track in SNI Area. 



28 




Figure 7. 



The San Nicolas Island Region. 




Figure 8. Mounting Arrangements 
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mast as shown in Figure 8. Thermo-Systems (TSI) Model 
1210 probes with sliding support shields were selected 
for these measurements. The sensor consisted of a plati- 
num coated tungsten wire 0,00015 inches in diameter and 
0,060 inches long, small enough to resolve the viscous 
dissipation scale without needing to apply wire length 
corrections. The wire was fitted with plated ends for 
isolating the sensing area and thereby minimizing flow 
disturbance . 

Electronics associated with each probe were a TSI 
Model 1054B Linearized Anemometer, a TSI Model 1056 Variable 
Decade Module and a Model 1057 Signal Conditioner. Schema- 
tics of the anemometer circuits appear in Appendix A. The 
anemometer had a linear frequency response from DC to 10 KHz 
and the variable decade module operated with a 0-60 ohm 
range. The signal conditioner had a linear frequency 
response from DC to 400 KHz and permitted voltage suppres- 
sion in one volt steps from 0-29 volts. The voltage output 
from the conditioner was recorded on magnetic tape. 

Sensor placement required 200 foot cables instead 
of the usual 15 foot cables. This resulted in a decrease 
in system frequency response but the decrease had little 
effect in the frequency band of interest. Ten probes were 
available for the experiment and they were placed into 
four groups. Each group of probes was calibrated with 
separate cables and electronics which permitted a complete 
system calibration for each of the four sensing levels. 
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The calibration was completed with a manometer and a TSI 
Model 1125 Calibrator shown in Figures 9 and 10. Empiri- 
cal relations used in interpreting the hot wire signals 
are presented in Appendix B along with the computer code 
used to evaluate these measurements and resultant calibra- 
tion curves for each probe. 

The performance of the hot-wire system was excel- 
lent throughout the experiment. Only one probe failure 
was experienced despite continued probe exposure during 
the data gathering period. 

2 . Temperature Measurements 

Temperature fluctuation measurements were made at 
two levels on the forward mast of R/V ACANIA as shown in 
Figure 8. Electronics used were TSI Model 1044 Tempera- 
ture Modules. The unit power supply consisted of six volt 
batteries. It can be utilized to make both temperature 
and AT measurements. The temperature module is designed 
for measurements of temperature variations down to 0.001°C 
at frequencies up to 1 KHz. The schematics for the Model 
1044 Temperature Module are shown in Appendix C. 

Probe failure occurred at the beginning of the 
first data run and to permit temperature measurements to 
continue the probe unit was fitted with a 60 ohm platinum 
wire sensor. 

3 . Parallel Measurements 

Mean measurements of the wind speed, temperature 
and relative humidity have been described by Cavanaugh 
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Figure 9. Probe Calibrator 





Figure 10. Vertical Probe Mounting with TSI Calibrator. 



(1974). Measurements of ship's motion have been described 
by Welsh (1974) . 

C. SIGNAL PROCESSING AND RECORDING 
1 . Analog Recording 

All data were recorded simultaneously in analog 
form with a 14 channel Sangamo Model 3562 FM Recorder*. 

All recordings were made at a recording speed of 7 1/2 
inches per second. At this recording speed the recorder 
had a frequency range from DC to 25 KHz. 

A schematic diagram of the complete recording sys- 
tem is shown in Appendix D. The amplifier and suppressor 
combination in series with the sensor electronics and re- 
corder channels were necessary to improve the signal to 
noise ratio of the inputs. The hot-wire signals were high 
pass filtered at 2 Hz during the majority of the runs. 

Two data runs were made without the filter to permit low 
frequency correlation between the velocity signal and 
ship's motion. The results of this correlation are dis- 
cussed by Welsh (1974). The temperature signals were low 
pass filtered at 500 Hz. In addition to the signal re- 
cordings, a voice track was utilized to record observa- 
tions pertinent to the data runs. Constant monitoring of 
amplified and preamplified signals was performed with an 
oscilloscope to evaluate the signals and to insure that 
the signal would not exceed amplifier or recorder input 
limits . 
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Real Time Analysis 



Real time analysis was accomplished using a Spec- 
tral Dynamics Corporation (SD) Model SD301C Real Time 
Analyzer. This unit was used in conjunction with a Model 
SD309 Ensemble Averager, a Model 13116-2A X-Y Display 
Oscilloscope and an X-Y plotter to provide an instantaneous 
evaluation of the hot wire and temperature signals. The 
equipment allowed positive determination that measurements 
were sufficient to define the inertial subrange. A system 
operation diagram is provided in Appendix E along with 
several of the X-Y plots obtained from the real time 
analysis . 
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. IV 



DATA REDUCTION AND ANALYSIS 



A. PRELIMINARY ANALYSIS 

The analog tapes were played on a strip chart recorder 
to permit initial evaluations of the data. From this, 
data periods were chosen for subsequent reduction and 
analyses. This preliminary analysis disclosed a recording 
problem that had previously gone unnoticed. A periodic 
tape transport vibration had occurred due to the 7 1/2 ips 
gear drive and this obliterated the recorded signal. For- 
tunately/ there were sufficient time intervals between 
vibrations to permit 20 minute data records to be obtained. 
It was from these data sections that the periods listed in 
Table II were obtained. 

The difficulties encountered in making the temperature 
fluctuation measurements, as described in Section III, 
were evident in the strip chart recordings. The tempera- 
ture data either failed to show sufficient signal strength 
or displayed a DC drift which rendered them impossible to 
digitize. An example of the two cases ilj shown in Figure 
11 . 

The paucity of temperature data led to a focusing of 
attention on the velocity fluctuation measurements. The 
criteria for selection of subsequent data periods for in- 
clusion in the reduction and analysis phase consisted of: 
(1) signal strength (2) averaging time (3) DC drift rate 
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Figure 11. Inadequate Temperature Signals Showing (A) 
Insufficient Signal Strength and (B) DC 
Drift. 



38 



and (4) availability of mean profiles. Sufficient signal 
levels were required to minimize amplification during the 
analog-to-digital conversion process. Twenty minute data 
runs proved sufficient to permit ample averaging time for 
Fast Fourier Transform ( FFT ) calculations. DC drift was 
kept at a minimum to permit ease in digitization and to 
allow greater confidence in the data samples. Finally, 
data periods were selected to correspond to times when mean 
wind information was available. A typical signal meeting 
the above mentioned criteria is shown in Figure 12. 

B. ANALOG- TO- DIG I TAL CONVERSION 

Analog- to-digital conversion was performed on a hybrid 
computer system located on the fifth floor of Spanagel 
Hall at NPS • The hybrid system consists of an analog com- 
puter, COMCOR Ci 5000, interfaced electrically with a 

> 

digital computer, XDS 9300. A detailed discussion of the 
operation of this computer facility was given by McKendrick 
(1972) which should be referred to for additional informa- 
tion. A block diagram, provided in Figure 13, illustrates 
the analog- to-digital conversion sequence followed in this 
s tudy . 

C. SEVEN-TRACK TO NINE-TRACK CONVERSION 

The analog- to-digi tal conversion process generated 
seven-track octal base data samples written on digital tape. 
The computer used to analyze the digital data was an IBM 



39 



3E 



w 



m 



S 

4 






ZsL 



Figure 12. Typical Velocity Signal. 
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Figure 13. A-D Conversion Procedures. 
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360/67 computer. This computer utilizes nine-track hexi- 
decimal based digital samples so a conversion between the 
seven-track and nine-track systems was necessary. A dis- 
cussion of the procedures necessary for making the proper 
conversion to the nine-track hexidecimal based system are 
described in Naval Postgraduate School Computer Center 
Technical Note #0211-08. The two primary programs involved 
in this conversion are the TAPEOUT and CONVERT programs. 
These programs along with a description of digital tape 
operations are discussed in Appendix F. 

D. DIGITAL ANALYSES PROGRAMS 

* Fast Fourier Transform ( FTOR) and spectral analysis 
(SCOR) programs were available in the IBM 360/67 subroutine 
library. The following description of the FFT prograrh 
package was taken from McKendrick. The programs were 
written by J. F. Garett and J. R. Wilson while students 
at the Institute of Oceanography at the University of 
British Columbia. The "OCEAN” subroutines of the UBCFTOR 
program instructs the IBM 360 to read the input tape, one 
block at a time, compute the FFT of the data block and 
then store the resulting Fourier coefficients on the out- 
put tape. The computer printout for the FTOR program is 
a listing of each digital channel and the Fourier coeffi- 
cients produced. An example of this output is shown in 
Figure 14 . 

With these coefficients, the UBCSCOR program is used 
to compute the spectrum for a given channel, or the cross 
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FTOR Output 



spectrum between a pair of channels, depending on what is 
selected by the user. The output of the SCOR program con- 
sists of a tabular computer printout, as shown in Figure 
15 f and a graphical plot of the spectrum versus frequency 
as shown in Figure 16. Appendix G of this thesis contains 
the preface comments of the FTOR and SCOR programs along 
with additional information for more efficient utilization 
of the spectral analysis package. Several difficulties 
were encountered during utilization of the analyses- programs 
and consideration of these difficulties is the primary 
reason for inclusion of Appendix G. 
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Figure 15. SCOR Spectral Statistics 
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V 



RESULTS 



A. ONE POINT SPECTRAL METHOD 

The procedure first selected for estimating £ , and 

hence u^, was to use single point values from the spectra, 
similar to that in Figure 16. A -5/3 slope was approxi- 
mated for the spectrum and the frequency and spectral den- 
sity, for a spectral value coincident with the line-, were 
taken from computer output similar to that illustrated in 
Figure 15. Equation 10 was then used to obtain £ from 
values of frequency, spectral density, wied speed and probe 
height. The computed value of £ was then used in Equation 
7 to obtain a value for the friction velocity (u*). The 
computer program used for these calculations appears in 
Appendix H. 

values were obtained by this method for all data 
periods and tabulated results for one minute and 10 minute 
averages are listed in Tables III-VII. This method ap- 
peared successful for spectra such as those shown in Figure 
16, which had a definite -5/3 slope for a majority of 
points. However, the: single point spectral method is quite 
subjective when used with spectra such as the one shown in 
Figure 17. It is easily seen that selections of unrepre- 
sentative spectral estimates could result in values of £ 
which would differ significantly. Because of the subjec- 
tive manner in which -5/3 slope lines could be drawn and 
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TABLE III 



Momentum Flux Calculations for 1103-1138 20 Sept 1973 
using the One Point Spectral Method 



Height 

(m) 


Time Average 
(min) 


e 

2 -3 

(cm sec ) 


u * 

(cm sec ) 


14.41 


i 


1.078 


8.727 




10 


1.319 


9.1.29 


13.20 


1 


1.010 


8.532 




10 


1.389 


9.109 


8.39 


1 


2.973 


9.871 




10 


2.075 


8.865 


4.52 


1 


2.418 


7.602 




10 


2.636 


7.812 
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TABLE IV 



Momentum Flux Calculations for 2008-2018 20 Sept 1973 
using the One Point Spectral Method 



Height 

(m) 


Time Average 
(min) 


e 

/ 2 ~ 3 
(cm sec 


u * 

) (cm sec ) 


14.41 


1 


9.742 


17.674 




10 


8.937 


17.272 


8.39 


1 


16.211 


18.291 




10 


11.679 


15.769 


4 . 52 


1 


20.806 


15.643 




10 


16 .128 


14.287 
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TABLE V 



Momentum Flux Calculations for 2023-2053 20 Sept 1973 
using the One Point Spectral Method 



Height 

(m) 


Time Average 
(min ) 


e 

2 -3 

(cm sec 


u * 

) (cm sec ) 


14.41 


i 


9 o 3 04 


17.505 




10 


9.556 


17.662 


8.39 


1 


14.042 


16.768 




10 


14.519 


16.956 


4.52 


1 


21.934 


15.829 




10 


22.634 


15.996 
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TABLE VI 



Momentum Flux Calculations for 2323-2341 20 Sept 1973 

using the One Point Spectral Method 



Height 

(m) 


Time Average 
(min ) 


e 

, 2 -3 

(cm sec J 


u * 

(cm sec ) 


•— 1 

• 

i— i 


1 


4.069 


14.925 




10 


3 . 543 


12.688 ' 


13.20 


1 


3.316 


12.973 




10 


2 . 576 


11.082 


8.39 


1 


7.965 


13.859 




10 


6.265 


12.813 
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TABLE VII 



Momentum Flux Calculations for 0637-0652 21 Sept 1973 
using the One Point Spectral Method 



Height 

(m) 


Time Average 
(min ) 


e 

, 2 -3 

(cm sec ) 


u * 

(cm sec ) 


14.41 


i 


1.923 


10.146 




10 


1 .114 


8.625- 


8.39 


1 


7.830 


13 . 991 




10 


5.660 


12.386 
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Figure 17. Velocity Spectrum Showing Deviation from 
-5/3 Slope. 
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a coincident spectral value selected, there was some 
doubt expressed in results from other than ideal spectra. 

B. INTERCEPT METHOD 

To avoid the subjectiveness of estimates of £ based 
on single spectral estimates, a second method was utilized. 
This method is based on Equation 10 rewritten in the fol- 
lowing form 

Ln k<|>(k) = -2/3 Ln k + Ln 0t£ 2//3 '(14) 

For this expression a plot of Ln kc})(k) versus Ln k would 

result in a curve with a slope of one and an intercept 

2/3 

value equal to Ln a£ . From this intercept, the value 
of £ can be obtained from 

Exp b 3/2 = 

L a J 

where b = intercept value. 

The advantage in using the i 
the single point spectral method 
determined from an average of th 
fined by the -2/3 slope. Figure 
curves obtained from the applica 
spectral results. Difficulty in 
was encountered in a few cases s 
in Figure 19 but in most instanc 
slope permitted easier detection 



ntercept method versus 
is that the intercept is 
e points along a line de- 
18 is an example of most 
tion of Equation 14 to 
approximating the slope 
uch as that illustrated 
es the forcing of a +1 
of spectral harmonics 
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Figure 18. Ten Minute Intercept Plot for 1103-1138 20 Sept 73 
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Figure 19. Three Minute Intercept Plot for 2023-2053 20 Sept 73 



which did not fit the -5/3 crit 
indicate a one minute spectral 
for that time average were simi 
clear that fitting a +1 slope t 
points would be difficult. 

The intercept method was us 
for all five data periods. The 
tions are listed in Tables VIII 



er ia . 


Table XI does no 


averag 


e since the 


curve 


la r to 


Figure 20. 


It i 


o thi s 


scattering 


of 



ed to calculate £ and u* 
results of these calcul 
-XII . 



a- 



C. COMPARISON OF METHODS 

A comparison of Tables III-VII and Tables VIII-XII 
indicates that somewhat higher u* values were obtained 
with the intercept method. The apparent validity of the 
higher values will be examined later. 

Considering the intercept method, a majority of cases 
yielded larger u* values for longer time averages. The 
importance of resultant higher u* values will be apparent 
when comparisons are made between shipboard results and 
results from other studies. The subjective nature of 
selecting spectral estimates in the one point spectral 
method is manifested by larger deviations in u* values 
associated with periods when spectra were not "near-per- 
fect." 

The one point spectral method was compared with the 
intercept method because of interests in designing pro- 
cedures for quick estimates of u* during measurements. 

If valid values of £ could be obtained from single spectral 
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TABLE VIII 



Momentum Flux Calculations for 1103-1138 20 Sept 1973 
using the Zero Intercept Method 



Height 

(in) 


Time Average / 
(min) 


e 

/ 2 -3, 

(cm sec ) 


u * 

(cm sec ■*") 


14.41 


i 


1.329 


9.151 




4 


1.378 


9.330 




10 


1.235 


9 .012 




20 


1.372 


9 . 248 


13.20 


1 


1.445 


9.139 




4 


1.379 


8.998 




20 


1.503 


9.260 


8.39 


1 


2.635 


9 .600 




4 


2.337 


9.224 




20 


2.622 


9 . 584 


4.52 


1 


3 . 086 


8.233 




4 


3 .148 


8 . 288 




20 


3 .689 


8.738 
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TABLE IX 



Momentum Flux Calculations for 2008-2018 20 Sept 1973 
using the Zero Intercept Method 



Height 

(m) 


Time Average 
(min ) 


e 

2 -3 

(cm sec ) 


u * 

(era sec ) 


14.41 


i 


10.288 


18.012 




3 


8.285 


17 .135 




10 


• 10.935 


18.474 


8.39 


1 


13.226 


16.437 




3 


11.828 


15.988 




10 


13.678 


16.622 


4 .52 


1 


17.765 


14.731 




3 


14.677 


13.846 




10 


18.127 


14.855 
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TABLE X 



Momentum Flux Calculations for 2023-2053 20 Sept 1973 
using the Zero Intercept Method 



Height 

(m) 


Time Average 
(min) 


e 

' 2 -3 

(cm sec 


u * 

) (cm sec ) 


14.41 


1 


12.071 


19.093 




3 


11.825 


18.985 . 




20 


11.540 


18.809 


8.39 


1 


16.723 


17.774 




3 


16.496 


17.694 




20 


16 . 243 


17.603 


4 . 52 


1 


25.023 


16.541 




3 


24 c 492 


16.432 




20 


24 . 149 


16.346 
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TABLE XI 



Momentum Flux Calculations for 2323-2341 20 Sept 1973 
using the Zero Intercept Method 



Height 

(m) 


Time Average 
(min) 


e 

2 - 3 

(cm sec ) 


u * 

(cm sec ) 


14.41 


3 


12.925 


19.533 




10 


13 .564 


19.849 


13.20 


3 


13.471 


19.235 




10 


14.009 


19.488 


8.3.9 


3 


14.418 


16.917 




10 


15.025 


17.151 
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TABLE XII 



Momentum Flux Calculations for 0637-0652 21 Sept 1973 
using the Zero Intercept Method 



Height 

(m) 


Time Average 
(min) 


e 

2 -3 

(cm sec ) 


u * 

(cm sec ) 


14.41 


1 


1.209 


8.868 




4 


1.955 


10.407 




10 


1.570 


9.675 


8.39 


1 


5.316 


12.131 




4 


8.451 


14 . 158 




10 


7.643 


13.691 
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estimates then the structure function measured may be 
valid. The latter has been shown to be a poor approach 
on the basis of these results. 



D. EVALUATION OF THE VARIATION OF U* WITH HEIGHT 

A long accepted concept in turbulence theory has been 
the existence of a con st an t- s t r es s layer. This concept 
is often used to evaluate the representativeness of accumu- 
lated data (e.g., Stegen, et al) and will be applied to 
test the validity of assumptions and expressions used in 
the present study. 

Table I was compiled by Stegen, et al for conditions 
when the measured mean wind speed was 5.12 m sec 1 . This 
mean speed corresponds to the wind conditions for the 
values tabulated in Tables VII and XII. A comparison of 
results in Table I and Tables VII and XII reveals lower 
values of u* in the present study. As previously men- 
tioned, the constant, (a), in Equation 10 is empirically 
derived and several values have been suggested in the 
literature. A value of a = 0.60 was used in compiling 
the results in Table I while a value of 0.48, suggested 
by Pond, et al (1966), was used in the present study. 

U* results listed in Table VIII support a constant- 
flux assumption as well as the validity of Equation 7 for 
u* calculations, for neutral or near neutral conditions. 

The anomalously low value of u* at 4.52 meters could be 
the result of acceleration due to the ship's bow. It is 
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seen that as the spectral averaging time increases the 
value of the flux increases. This would be consistent 
with an interpretation that increased averaging time 
would reduce the observed acceleration effects. 

Tables IX-XI correspond to periods of stable strati- 
fications while Table XII reflects a period of unstable 
stratification. It is seen from Equation 12 that neglect- 
ing diabatic effects would introduce a bias into the data. 
A stability correction could be applied to adjust the flux 
values for the more general case of stratified flow. 



E. 


THE VARIATION 


OF £ WITH 


HEIGHT 










Figures 


21-25 


show the p 


lots of 


Ln 


£ versus Ln z 


for 


the 


periods 


listed 


in Tables 


VIII-X 


II . 


On the basis 


of 


Equ 


ation 8, 


a plot 


of Ln e v 


s Ln z 


shou 


Id yield a -1 





slope for neutral conditions. Results from Stegen, et al f 
Figure 3, were in agreement with this prediction. 

In Figures 22 and 23 the slopes are less than -1 while 
in Figure 24 the slope is greater than -1. An initial 
concern when these results became known was whether possibl 
attenuation of the sensors at lower levels was a factor. 
Such a concern is logical since the sensors near the sur- 
face would be likely to lose sensitivity more rapidly as 
a result of increased salt loading at lower levels. 

The data period which occurred last is the one which 
shows the greatest slope which indicates that sensor deteri 
oration at lower levels was not responsible for the 
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Figure 21. Measured Values of e vs Height for 1103-1138 
20 Sept 73. For Scaling Purposes a Factor of 
5.0 was Subtracted from the x-Values. 
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Figure 22. Measured Values of e vs Height for 2008-2018 
20 Sept 73. For Scaling Purposes a Factor of 
5.0 and 1.4 were Subtracted from the X and Y 
Values, Respectively. 
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Ln C 




Ln Z 



Figure 23. Measured Values of e vs Height for 2023-2053 
20 Sept 73. For Scaling Purposes a Factor of 
5.0 and 2.0 were Subtracted from the X and Y 
Values, Respectively. 
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Ln C 




Figure 24. Measured Values of £ vs Height for 0637-0652 
21 Sept 73. For Scaling Purposes a Factor of 
6.0 was Subtracted from the X Values. 
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Ln £ 




Figure 25. Measured Values of e vs Height for 2323-2341 
20 Sept 73. For Scaling Purposes a Factor of 
5.0 and 1.4 were Subtracted from the X and Y 
Values , Respectively. 
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variation 



Further support of this conclusion 



in slopes. 

is provided by using Equation 11 as an independent method 
of determining the dissipation. 

For the 2023-2053 20 Sept 73 data period, the spectra 
at all levels were similar to that shown in Figure 26, 
where scale separation between the pure inertial transfer 
and the viscous dissipation regions was dl is cer nable . The 
upper limits of the inertial subrange were determined from 
the spectra and Equation 11 was used to obtain an estimate 
of the dissipation. In Table XIII the values obtained by 
this method are compared to those obtained from spectral 
estimates. The inner scale computation is independent of 
sensor calibration since it is determined by frequency dis- 
tributions and not on the spectral value. The apparent 
agreement between the two estimates is ara in situ calibra- 
tion which indicates the loss of sensitivity by individual 
sensors was not the cause for the deviation from the pre- 
dicted slope. 

The deviation from the predicted slope was also ex- 
amined for possible influence of non-neutral stratifica- 
tions or wind-wave coupling. Deacon (1951), from examina- 
tions of numerous vertical wind profiles, proposed the 
following expression for the wind shear 



8u_ _ ,z -3 

9 z kz 'z 
o o 



where 3 is a function of stability. 



( 15 ) 
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Figure 26. Twenty Minute Velocity Spectrum for 2023-2053 
20 Sept 73. 
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TABLE XIII 



A Comparison of Dissipation Obtained by Spectra (E) and 
Dissipation Obtained from the Inner Scale (e Q ) 



Height 

(m) 


e 

/ 2 -3 

(cm sec ) 


(cm 2 sec 3 ) 


u * 

(cm sec ) 


14.41 


9.556 


9.516 


17.662 


8.39 


14.519 


14.304 


16.956 


4.52 


22.634 


22.492 


15.996 
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The variation of 3 as a function or Ri , proposed by Dea- 
con, is shown in Figure 27. In Deacon's results, the 
variation of 3 may be described in terms of hydrostatic 
stability. For neutral conditions, 3 = 1 so Equation 15 
then reduces to Equation 6 and the -1 slope on the Ln € 
versus Ln z plot is predicted. 

For unstable conditions 3 is greater than 1 and for 
stable conditions 3 is less than 1. Deacon's results 
predict that the expected slope of a Ln £ versus Ln z 
plot would deviate from the neutral case with slopes less 
than -1 for periods of stable stratifications and with 
slopes greater than -1 for unstable conditions. Examina- 
tions of Figures 22-24 indicate why such considerations 
are important in interpreting the shipboard results. 
Figures 22-23 correspond to periods of stable stratifica- 
tions with the 2023-2053 20 September data period being 
the most stable. In both figures, slopes less than -1 
occur, with the more stable of the two periods showing 
the greatest deviation from the -1 slope. The results 
are, therefore, consistent with arguments proposed by 
Deacon. Figure 24 corresponds to an unstable period and 
shows a slope greater than -1 as predicted by Equation 15 

The results were further examined with respect to pre 
dictions from the Mon in-Obuhkov similarity theory. The 
theory states that in the boundary layer there exists a 
characteristic length determined by the stability and tur 
bulence which can be used as a scaling parameter. For 




Figure 27. Variation of $ with Ri (after Haltiner and 
Martin) . 
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> 






neutral conditions the appropriate characteristic length 

is the height above the marine boundary. 

Equations 12 and 13 show the effect of non-neutral 

considerations in developing an expression describing the 

vertical variation of e. Equation 12 indicates a relation 

ship between the scaling parameters Ri a rued L. Such a 

relationship was described by Busch (1573) in a review of 

information available in the literature. As presented by 

z z 

Busch, the function cj> (— ) - — in Equation 12 has the form 

ILL 

2 

(<J> - — ) where several flux-profile relationships have 

m L 

been suggested for 6 . Table XIV includes the formula- 

” m 

tions which Busch used to describe <J> along with constants 

m 

estimated by Paulson (1970). Re lat ion ships between the 

Monin-Obuhkov length and the Richardson number were first 

used to obtain stability length (L) corrections for the 

shipboard results, since Richardson numbers were available 

from profile measurements described by Cavanaugh (1974). 

2 

Values of L were then used to define (d> - — ) from reta- 

in L 

tionships in Table XIV. 

Figures 28 and 29 show the effects of stability correc 
tions for the periods with stable stratifications listed 
in Tables IX and X, respectively. Comparisons of the cor- 
rected and uncorrected results indicate that adjustments 
were toward the predicted -1 slope. The fact that the 
stability correction produces a -1 slope supports the basi 
premise of the Monin-Obuhkov theory and the validity of 
the relationships in Table XIV. 
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TABLE XIV 



Flux-Profile Relationships Describing the Stability 

Function, (f> . 

T m 



<f> 

T m 


Stability Condition 






(1-16|)" 1/4 




(l-12Ri) ~ 1/4 


Ri < -0.1 


(1- 3Ri ) _1 


-0.1 £ Ri £ 0 e 0 3 6 


0.88 (l-6Ri) _1 


0.036 £ Ri £ 0.1 
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Figure 28. Effect of Stability Correction for 2008-2018 
20 Sept 73. 
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Figure 29. Effect of Stability Correction for 2023-2053 
20 Sept 73. 
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The Monin-Obuhko v theory and the flux-profile relation- 



ships were also examined for the period with unstable stra- 
tification, Figure 30 shows the results of this examination. 
For -0.05 Ri -0.15 there is a tendency for adjustment 
toward the predicted -1 slope. However,, further decreases 
in the Richardson number caused an unexpected reversal of 
the adjustment. 

It is seen from Figure 30 that a value of -1.0 for the 

Richardson number results in a vertical variation of e which 

deviates significantly from the predicted slope. The reason 

for this deviation can be explained by examining the func- 
z z 

tion, d> Figure 31, from Garratt (1972 ), shows a 

m L L 

plot of this function using the expressions in Table XIV. 

To provide the correction necessary to obtain a -1 slope, 

z z 

((j) - — ) must decrease for decreasing values of — . It is 

m L L 

apparent from Figure 31 that this criterion is not met. On 

the basis of these results there is some evidence that the 

2 “ 1/4 

expression <p (1—1 6—) cannot provide the adjustment 

“ m L 

necessary to correct the plot for conditions of large 
instability . 

The fact that the stability corrections did not com- 
pletely restore the -1 slope, even in the stable case, 
leads to examinations of the possible effects of wind-wave 
coupling. Davidson and Frank (1971) observed wave-related 
velocity fluctuations and wave-related momentum transfer 
in spectral results obtained over natural waves. Davidson 
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Figure 30. Effect of Stability Correction for 0637-0652 
21 Sept 73. 
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Figure 31. Non-dimensional Dissipation Rate as a 
Function of z/L (after Garratt). 
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(1974) described the influence of wind-wave coupling as 
having the same effect as would be due to the influence 
of stable stratification. The wave-related effect de- 
scribed by Davidson was a reduction in momentum transfer 
from that associated with neutral conditions. The re- 
sults of the present study are consistent with Davidson's 
findings in that there is a decreased momentum transfer, 
u^, as the surface is approached. If this were due to the 
influence of wind-wave coupling, it would yield the ob- 
served deviation from the neutral slope, shown in Figures 
22 and 23. 



F. VARIATION OF TEMPERATURE FLUCTUATIONS 

Figure 32 represents the only temperature spectrum 
obtained from the five data periods which displayed the 
predicted -5/3 slope. Figure 33 is an example of the 
temperature spectra obtained from the majority of data 
runs. The paucity of temperature data precluded serious 
evaluations of the temperature fluctuation expressions. 
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L0G1O SPECTRUM (DEG/SEC /HERTZ) 




T J LEUEL 4 



Figure 32. Ten Minute Temperature Spectrum for 2323-2341 
20 Sept 73. 
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Figure 33. Ten Minute Temperature Spectrum for 2008-2018 
20 Sept 73 . 
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VI 



CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The results obtained from the one poirnit spectral 
method are generally lower than those obtained by the 
intercept method. However, the simplicity of the one 
point spectral method would permit a more rapid deter- 
mination of the dissipation and if a general estimate 
were all that was desired, this method would be accept- 
able. The results of the intercept method show excel- 
lent agreement among levels and with results from other 
studies. The longer time averages (-20 min.) appear 
to be necessary to avoid acceleration effects at the 
lower level. 

The close agreement between dissipation values de- 
rived from the spectral estimates and the inner scale 
calculations confirm the spectral estimates and also the 
choice of the constant, a = 0.48. Results on the varia- 
tion of e with height also verified the formulations 
used in calculating the dissipation and flux. In addi- 
tion, the sensitivity of the results to the varying con- 
ditions of hydrostatic stability and, perhaps, wind-wave 
coupling support the predictions proposed by Deacon and 
Davidson. 
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B. RECOMMENDATIONS FOR FUTURE EXPERIMENTS 

It is strongly recommended that a real time spectrum 
analyzer be used in future observational experiments. In 
addition to the capability for instantaneous evaluation 
of spectral characteristics, the resulting spectral plots 
could be used for general estimates of the dissipation. 

If observational experiments are to be made on a more fre- 
quent basis, there is insufficient time to properly 
evaluate the data before the next experiment. The use of 
real time analysis techniques would preclude the need for 
digitization and tape manipulations and save valuable and 
costly computer time 0 
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APPENDIX 
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APPENDIX B 



A, CALIBRATION AND SCALING PROCEDURES 
1 • Velocity Calibration 

Calibration of the velocity sensors was accom- 
plished using a manometer and a TSI Model 1125 Calibrator, 
The calibration procedure was to record voltage output ver- 
sus manometer height over a certain range of velocities. 

The computer program included in this Appendix was utilized 
to generate the calibration curve for further use in estab- 
lishing the calibration factor. This procedure was followed 
for each of the velocity systems and resulted in the cali- 
bration curves included in this Appendix. 

The calibration curve obtained from this program 
resulted in a plot of V versus /u where V is voltage and 
u is wind speed. 
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This curve describes a line 



v 2 = a(u) 1/2 + b (1) 

For use in data analysis it was necessary to obtain a 
scale such that 



C x e ' = u ' 



where C = calibration factor; e' = voltage fluctuation 
and u' = velocity fluctuation. Differentiation of Equa- 
tion (1) yields 



or 



or 



- 1/2 

2 V dv = 1/2 au ' du 



r 2x2 V u 1/2 , , 

[ ] dv = du 



4Vu 1/2 

[ ] e = u ' 



The calibration factor is then 



C 



4Vu 



1/2 



] 



where "C" scales between volts and cm/sec. 

The procedures for actual calculation of this 
factor consisted of (1) obtaining the slope of the cali- 
bration curve (2) selecting a "u" valve equal to that 
observed during the measurement period and (3) entering 
the curve with this "u" value to obtain the corresponding 
voltage . 
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2 • Scaling Procedures for Estimation of Epsilon 

In terms of the output of the spectral analyses 

programs the units associated with f4>(f) or kef) ( k ) were 
2 

(volts) . An examination of Equation (1) # Section II, 

2 2 

shows that the units should be cm /sec . Recalling that 
the calibration factor previously discussed scales between 
volts and cm/sec, it is necessary to multiply the E value 
by 






2 2 
cm /sec ) 

(volts ) 2 



-,3/2 



[ cm/sec ^ 
volts 



=' (C) 



2 2 

to obtain the value of e in cm /sec . 

To obtain the correct value of E it is also neces- 
sary to consider the gains used during the record and 
analog- to-digital conversion phases. In addition a scale 
factor is needed that scales between the analog voltage 
and the digital output. A 10 db gain was used when the 
raw signals were recorded and an additional 50 db gain was 
employed during the A-D process. The scaling between the 
voltage input and the corresponding digital output is ob- 
tained by considering that the maximum voltage input to 

the COMCOR Ci 5000 is 100 volts and this corresponds to 

23 

a digital maximum of 2 on the XDX 9300* 

The total gain and scaling for the observational 
data then becomes 



V (volts ) 



x G1 x G2 x 



SC ( 



digital 

volts 



Value obtained 
from IBM 360 



92 



where V - 

Gl = 
G 2 = 
SC = 

To obtain the 
be used with 



signal from velocity unit 

amplification prior to record (10 db) 

amplification during A-D conversion (50 db) 

23 

analog- to-digital scaling factor (2 /100) 

value which corresponds to that which can 
the calibration factor previously obtained 



V = (IBM 360 value x 1/SC x 1/G1 x 1/G2) 



The final scaling then becomes 



I 



[Scale x Gainl x Gain2] x C 



] 



or in this instance 

.23 



£[100/2“'' x 0.1 x 0.02] x C 



] 



C** #$*# $#*## 3{s ####$$ £:*.# £$*** $#*£# ## ^ ^c 4: ^ ^ 

aj<^ ^ £ * * # # ^ ^ ❖ ** 



c*** ** 

Q^C^C^t jjcjjt 

C*** THIS PROGRAM IS DESIGNED TO COMPUTE A CALIBRATION ** 
C*** CONSTANT FOR USE IN THE SPECTRAL ANALYSES PACKAGE. ** 
C*** THE INPUT PARAMETERS CONSIST OF THE MANOMETER HEIGHT** 
C*** AND THE CORRESPONDING VOLTAGE. THE PROGRAM WILL ** 
C*** CALCULATE AND PRODUCE A PRINTER PLOT OF VOLTAGE ** 

C*** VERSUS VELOCITY. A BEST FIT CURVE CAN THEN BE ** 

C*** DRAWN THROUGH THE DATA POINTS TO OBTAIN THE ** 

C*** CALIBRATION FACTOR. ** 

C*** ** 

Q^c^c jjcjjc 

C*** STATEMENT OF VARIABLES ** 

C*** ** 

C*** H MANOMETER HEIGHT IN INCHES ** 

C :&;& Ox. *ic 

C*** V VOLTAGE ** 



C*** 

C*** 

c * * * 
c*** 
c*** 
c*** 
c*** 
c*** 
c*** 
c*** 
c*** 

c ** ** ** ^^**+-* £ 

0** 'V' 'V' 'I' -v- '1' 't' 'l' ~l~ -p 'i' -V- 

o*************** 



VSQRD. 

U ..... 

X 

USQRT . 



VOLTAGE SQUARED (ORDINATE OF GRAPH) ** 

VELOCITY IN FEET PER SECOND ** 

SQUARE RCGT OF H ** 

:};5jc 

SQUARE ROOT OF UtABCISSA GF GRAPG) ** 

U# *JU JL X w %L> JU <J U Jy Vf ■JU vL JU ^ «JU< JL Jf >1^ U# JL <JU ^ JL JL -JU JL JU JL> X >u <JU -JU ^ X X X X 

x ^ ^ x *^r T 1 ^ -y* ^r* ^ -v* -v "Y' -v* ^ -y* x* if 1 •p “*p x ^ ^ ^ x ■*Y k 




J. JO U* jL. JJ -l . X JU 



t* -r* -nr -r* -r> - * v n' r T 1 o'T 1 -r m -r -i* vr *v »r» -r* *r 



REAL*8 TITLE( 10) 

DIMENSION VS QRD ( 15 ) , JSQRT ( 15) , U ( 15) 

READ (5,20) TITLE 
20 FORMAT (I0A8) 

WRITE (6,30) 

30 FORMAT ( • 1 1 , T9, ' H 1 , T1 8 » ' V* ,T28, ' VSQRD' ,T39, ' U' ,T47 , 'US 
1QRTM 
N = 15 

DO 60 1=1, N 
READ (5,40) H , V 
40 FORMAT (2F10.3) 

X=SQRT ( H } 

VSQRD ( I ) =V**2 

U(I )=65.3*X 

USQRT ( I ) = S QRT ( U ( I ) ) 

WRITE <6,50) H,V.VSQRD(I)»U(I ),USGRT(I) 

50 FORMAT ( • O' , 5F10.3 ) 

60 CONTINUE 

WRITE (6,70) 

70 FORMAT Cl') 

CALL PLCTP (USQRT, VSQRD, N, 0) 

WRITE (6,30) TITLE 
80 FORMAT ( • O' , 10A8) 

STOP 

END 
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Figure El. Real Time Spectral Analysis of Velocity Signal 
for 2008-2018 20 Sept 73. 
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Figure E2. Real Time Spectral Analysis of Velocity Signal 
for 2023-2057 20 Sept 73. 
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APPENDIX F 



A. TAPEOUT PROGRAM 

TAPEOUT is a general purpose program used to determine 
the label and file contents of magnetic tapes. The program 
can take an otherwise unreadable tape and provide the user 
with sufficient information to permit access to the tape. 
TAPEOUT will print a file number, density, parity, length 
in bytes of the record and a dump of the record. In this 
study the two key outputs were the record length and record 
dump. These outputs permitted the user tho obtain the block- 
size and ascertain the location of erroneous records. The 
following is an example of the TAPEOUT deck set-up used in 
this s tudy „ 

//JOHNSTON JOB ( 2 4 1 3 , 0 3 9 2 , XM2 4 ) , ' JG»HNSTON SMC 2561' 

//EXEC TAPEOUT, PARM= '2.0, , 2 1 

//TAPEIN DD UN IT=2400-1 , LABEL= ( , BLP )) , V0L = SER=KARLYN 
/* 

TAPEOUT is part of the IBM 360/67 subroutine library 
and a discussion of the parameters (PARM= N ) is available in 
the computer center consultant's office and also in Com- 
puter Center Technical Note No. 0211-08. 

B. CONVERT PROGRAM 

The seven-track octal tape is converted into nine-track 
hexidecimal data samples using the CONVERT program. The 
block characteristics of the data samples are still 
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maintained on the nine-track tape but two additional para- 
meters are affixed to each block. These two parameters, 
NCHAN and KMAX , are affixed to the beginning of each block 
NCHAN is the number of channels of analog data digitized 
on the seven-track tape. KMAX is a numerical value set 
equal to the maximum number of samples per block. The in- 
sertion of these numerical parameters at the beginning of 
each block and the change in the block format by the FORM 
subroutine will put the seven-track tape in the proper for 
mat for use in the UBCFTOR program. The CONVERT program 
used in this study along with sample JCL cards is given 
below . 



//JOHNSTON JOB ( 2413 , 0392 , XM24 ),' JOHNSTON SMC 2561' 
//CONVERT EXEC FORTCLG , REGION . GO=100K 
//FORT. SYSIN DD 

DIMENSION IDAT ( 1024 ) , DAT (1024) 

FACTOR = 1. 

REWIND 2 
REWIND 4 
NCHAN = 8 
NRECL = 1024 
J = 0 

10 READ (2, 20,END=60,ERR=40) IDAT 
20 FORMAT (20 (50A4) ,24A4) 

J=J+1 

CALL FORM (IDAT, NRECL) 

DO 30 1=1, NRECL 

WRITE (4) NRECL, NCHAN, DAT 

GO TO 10 

40 WRITE (6,50) J 

50 FORMAT ( ' 0 ' , 5X , ' READ ERROR, RECORD N0.=',I4) 

GO TO 10 

60 WRITE (6,70) J 

70 FORMAT ('0',5X,'END OF TAPE , RECORD N0.=',I4) 
END FILE 4 
RETURN 
END 

/* 
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//GO .FT02F001 DD UNIT=2400-1, VOL=SER=KARLYN ,LABEL=(,NL) , 
// DISP = OLD, DC B = ( DEN= 2 , RECFM=VS , BLKS I ZE= 8 1 9 2 ) 

//GO.FT04F001 DD UNIT=2400, VOL= SER=NP S 1 7 0 , DSNAME= DEREK , 
LABEL= ( 1 , SL ) ,DISP= (NEW, KEEP) , DCB= ( DEN = 2 , RECFM= VS , 

// BLKSIZE=8200) 



It is noted that the user 
of the dimension of the arrays 
of NRECL • IDAT and DAT should 
x N CHAN thus signifying the ne 
ing parameters. The argument 
value equal to NSAMP x NCHAN • 
tor has been set equal to one 
found to be easier to perform 
in one step following the use 
grams. Computer Center Techni 
additional information concern 



inputs to the program consist 
IDAT and DAT and the value 
be dimensioned at least NSAMP 
ed for recording the digitiz- 
NRECL must be assigned a 
The value' of the scaling fac- 
in this study since it was 
all calibrations and scalings 
of the spectral analysis pro- 
cal Note No. 0211-08 provides 
ing the CONVERT program. 



C. DIGITAL TAPE OPERATIONS 

Unlike the Hybrid computer system, where all tape mani- 
pulations were performed by the user with technical assist- 
ance from the computer facility staff, the IBM 360 computer 
system manipulations were performed by the staff members. 

In order to have the tapes handled as desired, explicit in- 
structions had to be given. These instructions were given 
in two ways. The first instruction consisted of a Job Re- 
quest Form addressed to the computer center staff identify- 
ing the tapes to be utilized and whether they were to be 
read off of (RINGOUT) or written on (RINGIN). The second 
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instructional procedure was conducted usimg job control 
language (JCL) addressed to the computer and read as part 
of the submitted program. 

1 . Job Control Language 

Job control language cards providce information to 
the computer about tape identification, tape disposition, 
tape data files, tape mount numbers and data format. The 
cards are appended to the program and rea<d through the "hot 
card reader" as part of the program, whereas the Job Request 
Form is submitted "over the counter" to tdhe computer center 
personnel. An example of a JCL card for itape processing 
would be : 

//GO.FT08F001 DD UN I T= 2 4 0 0- 1 , V0L=SE R=KARLYN , LABEL= ( 1 , SL ) , 

// DSNAME = DEREK, DISP= (NEW, KEEP ) , DC®= ( DEN= 2 , RECFM=VS , 

// BLKS IZE=8192 ) 

While the majority of the JCL card remairas unchanged, there 
are portions which are subject to change Ifrom tape to tape. 
These portions are discussed below. 

FT08 - This group indicates the tape mount on which the 
tape is to be placed. The two digit number (08 in this case) 
must match the number provided in the program's "READ" or 
"WRITE" statement. 

F001 - This group indicates the sequential number of 
passes through the tape which has occurred up to this point. 
In this case it is the first pass. 

UNIT=2 4 00-1 - This group specifies the tape as a seven- 
track tape. The nine-track designation is 2400 o 
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VOL=SER=KARLYN - This group identifies the tape to 
be processed. In this case the name "KAtLYN" signifies 
the tape is an "external" tape with the name "KARLYN" 
physically appearing on the tape. An "internal" tape 
would have a serial number such as NPS2L1. These tapes 
are kept on file in the computer center and cannot be 
physically removed from the facility. 

LABEL=(1 / SL) - This group identifies the tape file to 
be processed and specifies that the tape has a standard 
label. "External" tapes may be processed as no label 
tapes (NL) but "internal" tapes will be given a standard 
label before the user gains access to them. This is accom- 
plished by the computer center staff processing a "canned" 
program which places a standard label on the serialized 
tape . 

DSNAME= DEREK - This group specifies the data set name 
(DSNAME) applicable to the specified LABEL group. The 
data set name may be a unique name limited to six alpha- 
numeric characters. Once a file has been given a data set 
name access to the data file can only be obtained if this 
same DSNAME is used. This process can be circumvented by 
the use of by-pass label processing ( BLP ) but does provide 
protection against accidental erasures. 

DISP= (NEW , KEEP ) - This group specifies the tape dis- 

position once it has been processed. 

DCB= (DEN =2 , RECFM = VS , BLKSIZE = 8192 ) - This group is re- 

ferred to as the data control block. The density ( DEN= 2 ) , 
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tape record format (RECFM=VS) 4 nd block size (BLKSIZE= 
8192) are stipulated using this group. The density speci- 
fies the number of bytes per inch, in this case the 
number 2 denotes 556 bytes per inch. This corresponds to 
a manual setting of the DENSITY dial on the face of the 
seven-track tape unit. RECFM=VS specifies that the record 
format is variable spaced, the result of recording and 
digitizing on the seven track unit* The block size speci- 
fies the number of bytes per block of data. If 2048 words 
per block were used, four bytes per word would result in 
8192 bytes per block. 
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APPENDIX G 



A. GENERAL 

All the UBC time series analysis programs ( FTOR , SCOR, 
FC PLOT) were originally listed on faculty tape NPS216, 
file 1-3 respectively. The programs have since been placed 
on computer center disk, MARY, and are stored in machine 
language which disallows a listing of the program for user 
study. However, source decks for UBCFTOR and UBCSCOR have 
been retained by Dr. Kenneth Davidson of the Meteorology 
Department and by Dr. Noel Boston of the Oceanography De- 
partment. These source decks can be used to obtain a 
listing of the programs for individual study. The most 
recent computer center documentation on the use of the time 
series analysis package is given in the following memoran- 
dum written by Miss Sharon Raney of the computer center 
staff. This memorandum provides the information and JCL 
necessary for utilization of the disk-stored analysis 
programs . 

B. FTOR 

The preamble to the FTOR program is given following 
this discussion. It provides the user with the informa- 
tion needed to prepare the input data cards necessary for 
program execution. One of the most misleading parameters 
called for in the input values is the variable NBLOCK . 
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This value is provided by the user on the second data 



card with the value supposedly entered in columns 4-5. 

In actuality the value is read as a five digit integer 
thereby making the maximum value of NBLOCK 99999 vice 
99 as implied by the program preface. 

An important consideration involves the value of the 
variable MPRINT . This variable controls the computer 
printout of the FTOR program. If the user desires to sup- 
press the printout the insertion of a specific value can 
save considerable paper and cut computer time significantly. 
In the unsuppressed mode, the output will follow the format 
shown in Figure 13, Section IV. The suppressed output will 
consist of one line outputs for each data block which will 

inform the user that "coefficients for block have been 

written on the file of the output tape." By using the 

suppressed output, computer turn-around time was signifi- 
cantly reduced. 



C. SCOR 

The preamble to the SCOR program also provides the user 
with the information necessary to compile the data input 
cards. The SCOR program requires the input of several 
variables whose values are determined during the A-D con- 
version process. One of the values on the first data card 
is the number of data blocks on which the spectral analysis 
will be performed. This value is set by the user and can- 
not exceed the number of data blocks processed by the FTOR 
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program. As was the case in the FTOR program the maximum 
number of blocks which can be handled is 99999 and not 99 
as the preface would lead the user to believe. 

As was discussed earlier in this thesis the output of 
the SCOR program consists of a computer printout and a 
graphic plot. The graphic plot consists of an output in 
which spectral density is plotted against frequency for 
all channels being processed. Since the original formula- 
tion of the .SCOR program, the computer facility at the 
Naval Postgraduate School has installed a new plotting 
routine called NEWPLOT. This routine limits the user to 
30 inches of plotting paper. The author found it necessary 
to submit a job request asking that plotting routine CAL- 
PLOT be used to prevent early termination of spectral 
plots consisting of four or more channels. The computer 
center staff has been alerted to this problem and any 
future user should consult with a staff member before 
using the SCOR program. 
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8 June 1973 



MEMORANDUM 

From: Sharon D. Raney 

To: Users of UBC Time Series Analysis Package 

Subj : New Operating Instructions for UBC Time Series Analysis Package 

1. The new FTOR deck setup is given in Figure 1. 

The dimension size of array DATA must correspond to the block size 
of the data on the input tape. This tape is generally the data as 
converted from the seven track XDS generated tape to the nine track 
IBM/360 tape. 

Formerly, this array was dimensioned in subroutine OCEAN 2 and had 
to be changed whenever the sample size changed. This meant the 
entire source deck for FTOR had to be used for each run. 

FTOR is now compiled in FORTRAN IV (H) which generates a more 
efficient code and consequently speeds up execution time. 

2. When the same dimension size is to be used repeatedly. Figure 2 
shows how to obtain an object deck and how to execute FTOR with the 
object deck. 

3. Figure 3 shows the deck setup for SCOR. 

4. Notice the change in the format of the Job Control Language (JCL) 
to define the tape files. These changes enable more efficient 
processing by the operator. Note the use of the IN parameter in 
the LABEL parameter. This must be specified when the file is to 
be read only. OUT is used when the file is to be written only. 

All data sets are PASSed not KEEP. The UNIT— AFF is used on sub- 
sequent files of a tape as is the VOL=REF. 

5. Any questions concerning these new operating instructions should 
be directed to Sharon Raney at X2631 or Ingersoll Hall Room 107. 
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FTCR DATA CARDS GO HERE 



//FTOR EXEC FORTHCD, REGION. F0RT=150K 
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Figure G2. 



115 



// DISP«( NEW, PASS) » LABEL* (2, SL, ,OUT) , DCB=( DEN = 2 , R ECFM=VS , BLKS I ZE= 12300) 
//GC. FTC3F003 00 UN IT = AFF= FT03F00 It V0L=REF=*. GO .FT03F001, DSN=C0EFQ1 , 

// D ISP = (NEW, PASS) , LABEL = (3,SL, , OUT) , DC8=(DEN = 2 ,R ECFM=VS , BLKS IZE = 12300 ) 
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CARD TERMINATES THE RUN 









-FTOR- 










> 


3 








CQ 


O 










2 








Q ^ < 


3 








3 O 


3 








3 0C Z 


Of 








030 3 


3 


OO 






312 


OO 


3C 






31- O 




O 






O h- 3 


< 


3 






LLOOOO O 




3 






30C O 


3 


o 






«-»3 *~ 1 } 


O 


LLO 






CMZ3 O 




LU 






200 


O 


003 






O i£<0 30 


2! 


<3 






3 0X3 cl 


3 


< 






3 OOZ 33 




00 O 






3 3 0 o:m 


3 








< 002:2: •< 


X 


00 






0 *— <<q:— 2 : 


f— 


QQ 






X<O00M 




LUQC. 






x uoxq:o< 


h- 


2:0 


O 




O 1— 031— 


< 


CCjZ 


LU 




< WSQZ 




0 


5T 


— * 


3 OUJm ^0 


3 


U-nO 


cC 




— IX 0 0 


a 


00 m 


O 


— < 


00 l-|-Z~ 


< 




3 


<<M 


0 a oz>- 


f— 


< 


ZL 00 


<w 


Q C Z 31 — O < < 




CL 3 


O Z 


O' 


0 G<3iz 


i— 


O 


O < 




12 0 L/)0 


3 




3 OL 




HX 


3 


*:> 


QQ »— 


— % 


vO ^UU.h'K' O 


1— 


o< 






LO 0<3QCf\|CL 


3 


Oa: 


00> 


#• 


CM 2: 33 <w 0 


O 


-Jo: 


00 JO 


-^r~4 


cl Oa + s 




C0< 


UJ^LU 2 




3 <a:o 3 


3 




3-^Z3 


wOO 


0 XO > 3: 


X 


XZ 


Q_OZ3 


2:|- 


LL/^friUjJ 


f— 


o< 


Z3<o 


<HH 


00 CL 3 <03 




< 


<3X3 


ZZ 


O LUOhZ 


2! 


UJOO 


alOOQQOO^ 


-13 


q: oza'M^h 


O 


*— « 


3 3 


*:o<— 


O 3 *— * 3 O 3 LO 




a: 


CDILULIU. 


w< NZ 


O 1 — 2 o<o< 


Z 


OQ 


5Z03s£)(J) 


z II w 


3LLZH3ZZJ 


3 


3a: 


3 00 Z 


< II -J 


aco*— • < »—»3 


h- 


O 


Z^wQlHZl ->< 


1 — 3 33 


1 — 


LUO 


a: 3 33zu-^u 


3hozxajax 


1 — 1 


CLLU 


3 ^ CO Cl CQ Cl < — — 1 


<00ZOI— X3H- 


3 


<tcc 


00 USK^SI 


r-H | II 


Ummj 1 — LO 




h- 


DODI-D<h 11 | ^ 


•-•OOZ 3 h- 






Q3ZZZ<^ 




OZmooO>LU< 


OO 


l-< 


*— CO || 


— wCNJ^Z 


00<Q cOh-X 


»— * 


30 


II II II II 


4- + 


— IOI — Odl — 3h 




Q-*-h 


II II — 


+ + 30 


ZOo^OJ 


3 


1-0 


• — ■ ^ ^OO^N't 


11C02<^G-G 


— loO 


30 


»h c\i ro < un 0 h —1 ro 3 ,-h 


< Zo Q.2Z t/) 


HHNZ 


03 




IW w W w 


h-< 2! OO 


30 




33333303303 


Z 0< 3 lO>* 


a 


U-J—LUCLCLClclClClClCLClClCl 


>3X<3 < 


33 


O00CL<C<<<<<<<<< 


1— <0030>a: 


03 


a: < 1 — | — | — 1 — | — 1 — | — 


1 — 1 — 1 — 1 — 


X^2Zq;3J(1 




h- CLCLCLcLcL cLCL<LiLCC CL 


33 -ft Z< 


33 


<3 cL <<<<<<<<<<< 


Z<— I'M Oh 


2!0 


< 






h- 


3 


cL 










0 








z 


LL 








< 



• • 

ilcl 

• • 



O OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 



GARRETT INSTITUTE OF OCEANOGRAPHY 

WILSON UNIVERSITY OF BRITISH COLUMBIA 

VANCOUVER 8, CANADA 



SCOR 



OO OOO O OOOOOOOOOOOOOOOOOOOOOOOO oooo oooooooooooooo 

— K\icn<r Lnsor^-cooo— <rm or^oo'0-*rgcn<u^oi^co^o-ir\jcn>tioor-.oo 

OnOOOOoOOr-f^HrH^r^ r -iHHH^N^^(\J'N^N(\J^^^r0^rr) mrOrOrnrO^ n 4 -^ -1" <r-4‘'T-t<r 
OOOOOO OOOOOOOOOOOOOOOOOOOOOOOO OOOO OOOOOOOOOOOOOO 

ooooo OOO OOOOOOOOOOOOOOOOOOOOOO OOOOOOOOOOOOOO oooo 







U2 


XO 














o 








-JO 


<t 2 


X 












2 Ol— 








Q <t — ♦ LLI 


CDC 


LU 












•— » 2uO 








ZXHUJ 

Qt< COO 


I •» 


h-CQ 

22 












1— •*-^X 
31 — 2 








O 2 — ^ >- 


CQ 


U3 






2 






X02ll.lL 








cQ 


l-LU 


h2 






C 






ZZ H o 








UX 11)0 • 




O'-" 






X 






o CDlU 








hOOUJ>2 


Ol-X 


— • h— 






CD 






OXU X LU 


Q 






>00 I/O b- LLI 


NO(- 


LL X 






O 






OCQh-U 


LU 






CQLL) 


• — * — J 


LLO 






X 






LU »— i 


X 






Ql>- CX 00— • 


X X LL 


LUOO 






xoo 






XQUoou 


►— * 






OOOCCoo2CD 


o o 


o 






c 






h-UX^ 


oo 






LU 2 O XU 


XU 


o\ 






oo 






Ol— x 


LU 






OoO 2 OQ 00 


h-2 


u 






K-.2 






X 3 oo 


a — 






DUQKO *-< 


c<o 


X 1 — 






XO 






OQ22C 


o 






oooh-uu 


Sm 


LU C 






h-— ‘ 






uoocu 


c^ 






03 oocbc-^ 


Jt-th 


•-<2 






1 — 






X 


i — 






xq>- occm 


-Jh-c 


X — 






X 






1— XQcOX 


ex 


(X 




CLOO ►XHw 


<00M 


3l- 






XO 






3 uoo 


QC 


u us: 




X2UCO 


u> 


OOO 


X 










OlUJUOIL 


2 


m CQ 


h- 


LU CL LU O 2 LU Q 


2 U 


LL UJ 


u 




QLL 






h-X3 


X w 


X»~< 


1— 


CL OCcdclz: 


k-*LUO 




oo 




LU> — i 






3CQU|— 


o 


3 


Hi 


<iut?a:^M 


X 


XU 


z> 




OOh- 






O 0X2 


UQ 


02 


x 


hOOULUi) 


KQ 


LU> 






3^ 






U)Xh- u 


UJ 


IL-O 


X 


LUX > LULL 


• X 


_)—• 


> 




u 






JhX oo 


XOO 




e 


2XlL<rXOH 


U2< 


XCD 


CQ 




U3 






COU2U 


LU 3 


2 : 5 : 


cd 


Oh h- w 




ST 






X^ 






oux»-*x 


cO 


o< 




X X2 UJ 


CO 2 


QO 


Q 




c 






•—•X X 


2U 


lXQC. 




LUU|-<X32 


coc 


Uh 


u 




Q 






X uu 


3 CQ 


LL CD 


X 


O^H LU_J'JJ 


jUJh 




1— 1 




tn 2 






a xcoooo 


2 


O 


o 


</) IS •>C LU 


hQOO 


UQ 


u 




h-C 






C LU 2 ►— • 


O 


OOX 


“D 


1 — oo >0>3 


C3 


Xu 


X 




►— » 






xxuo 


21— 


UQl 




zs:qu j— 


>JQ 


H-O 


X 




200 






CD|— u t/) « 


CJ 


»—4 




U3 U 2 2X U 


cou 


O < 


Z) 




3r— 






LUm>QO 


►— »t/) 


|— » 




hoOXJQu 131 


O Zh 


• t/)X 


oo 




2 






xx^uxu 


1 — — 1 


OOX 




O mZm< 


• OO 1 — 1 0 


(\)> — i u 






V-U 


oo 




D5 2 Coo 


c u 


_<0 




m'JJilGlO'IZNm LU 


> 


u 




3*— 1 


X 




o O O OOO 


02 


h- h“ 




Uh- LUOO Z5 


UX 


c 


00 




xo 


LU 




>HX u 


•—•2 


CU 




LXC LU X U CY X \ h- < X > 








H-h-. 


1 — 


Q 


2COOLLO 


U< 


h-- 




LU i — •CQU-Qf'OI — 


— «Z5 > LL! O >- 


— — oo 


h- 




3LL 


UJ 


LU 


>UOXQQ 


' X 


00 


cr> 


O X CD LL O • 


-^XX 2C0 


f\ J-'X 


OO 




OU 


2 


X 


UOO C X 


1—0 


>- 


vO 


OOL>3ijJ LU '\J 1 — UJ LU LU 


Wr-jO 


3 




XU 


C 


•— » 


• — ' LU X) LU I— Ol 


2 


2TC0 


O' 


OCOX2X 


-3hXD2 


— 1 o 


2 




h-o 


X 


3 


I — X»— « X LU 


UU 


3 




XX2-J LUc/) 1 — ‘021 — C? 3 


— I 






30 


c 


O 


2X1 XO0LU 


QO 


(XO 


•* 


LUX LL > O -ft UX 


# H-CO 


oo 




O. 


X 


u 


U 1— C CO 


• — ♦ 


1— LU 


h- 


hQ.QOO*-«U 


— O OOXH* 


— X 


UJ 




2CDX 


X 


QU CD 


X 


oo 


<\J 


X <1 LU LU CD 


ULU>-ULLO 


r-4 C?LL 


2 




►—•20 1 




h aQLL Q 20 


xu 


ULO 




3 OO— »h- U 


wh-OH* u 


WOOD 


t — i 




*-Hi 


HOH 


oo 


OXi-nh- 


U CO 


QlO 


> 


OLU3LL2UX 


HhZ< h-X 


X,ft 


1 — 




U>LLX3 


) — 1 


Ol- 


2002 


00O 


x 


XX ►— ‘LU Xh— 


O LU O C OO 


— , • 


3 


U 




H— X 




uOOOUOfO 


233 


Xh- 


c 


hiOO^< 


+ JDhI 


1 2> 


CD 


u 


OCL><- 1 — 


1 — 


32 Jh 


3 2 


OOCLvO 


3 


oo X Xlo 


ClOOH-J 


wLU 


OL 


00 


i — iq.cCjU3 


3 


QZZOO JZ 


— j II 


00 


2 


Q5DQlm-> 


— U2 X 


^h-»Q 


CD 


c 


CD3 


OO 


X 


LU j » — i O LU 


O II 


OULU 


<1 


<O(/)U0LLZONlLo:^h 3 


i — f * • 


3 


u 


Ooou 


2 


alhJXiLH 


ox 


X X LU 


“D 


UX LL 02 


wQLL <h“ 


W 


oo 


► 


u 


LU— » O 


• — » 


hhQI — o 


LU X 


o<o 




XUUClU— 'UJX X < 


— + 1 




o 


— * 


OOO LU 




3 O U*-< 


2oo< 


i_o 


2 


X f.l-3-ft >UC^X 


-ft OOOO 


CD 


» — 1 


OJJDQh 


CD 


C ? • • X LL 


—0 2 


O 2 


O 


2 01— 20X0 


-h-ii^za 


X 




h- 


2x0X2 


2 


LU LL C\J | — 1 — U 


►QO 


Zh 


» — ♦ 


C 2 O ^'iJHLUh W C 


r\) 2 II O 


HHOO ► 


» — i <i o c * — 1 


> — i 


x LLr — o u 


Qhm 


<zo 


OO 


X C • X rt X X 


W~ U Z) 


— o 


^uu 


3HXOX 




< u «o 


X 


L±J vO 


w 


CD <XILCOLLXXOO<00<3 


XX u 


Ooo< 


Ow 


X w x 


O 


ool— h— X CjO 


c IA 




> 


O »*X U(Z) 


C LUO U 


UJcC 


uc> 


u 




u 


*— 'OOt/) X 


O —* 


OO 


LU 


X 2 1 — OO 2 LU X 


«*— - ►— • m LU 





uxu 


U II 


M II 


u 


XQC<X 


0 s 1 


X — *2 


X 


X C O LLI > — iQ(J 


CJ-XX 


C LL 


oxx 


OfO 


IX o 


o 


Q LU ►— • 3 O U 


h- | <■ 


h-U-LU 




XUJ>U c 


c o X u 1 — 


UJ20 


LL 




LL 




LL 


XalUO *-• 


OO—i — 1 


O LL 1 — 


h— 


OOCDXC Xu 


— lu n 


X — 1 












Cl— >-< X 


X 


LLI LL! L/) 


OO 


MQoOhQO 


Q> Q 


U 1 — 


U 




U 




u 


02U 1 — 3 


> — ( 


QlO> 


C 


X X O O 2 1 — 1 — 


2 2 


X oo 


X 




X 




X 


lillu.<u 


u 


</iO 00 


u 


1— XOOCLLC 


•—e c 


I£U 


h— 




1— 




1— 


C Ol — O QU 





121 



-SCOR- 



OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
Or-icMrovr«A vor^co o^o^(Mro<r:nvcr- -oo'Or-tfNiVi ;rovOHooo.-injrns^-in<)Hcoo>o r - H c\jr<v;f*inv 0 
inin^iiriLnLnin^iriinvOOvOOO 0^000 DMHHMHH^HHMaooooocorocxicocoac>aQa' 0 ''a>;j><J'vJ'a' 
OOOOO 000 o 0000000 00 000 000 000 OOOO 00 00000000 0000 00 
00 OOO OOO O OOO o 0000 00 00 OOo 0000 000 00 0000000 0000 000 



0 


> 

O 






















O' 


2 






















sr 


UJ 








0 














O UJ 


DO 






0 


GO 




, — , 










O CL 


Oo 






< 


O 




• 


X 




—a 




< 


>LUO 






UJ 


+ + 




UJ O 


O 


• 


Hoi 




' H 


OCLQ 




0 


H 


O 




CD < 


2 


Of 


2UL 






2U 1 




LU 


00 


OO 




3 


0»~* 


OJ 


LL —* cC 




H 


UJ 




H 


2: 


O 




_J O 


O 




OOl — 


•» 


2 


=)hh 




H 


*— • 


LL 


• • 


_JLU 


<QC 


OLL 


CL oO 


< LU 


LU 


CfroOoO 




O 




ai 


00 


m CQ O 


UJ 


—I 


H 


CC OC 




lu22 


O 


O 


O 


0 


—lO 


2 2 


HCL 


LUO 


2 — 1 LL 


H C 


o< 


O^mhh 


LU 


CL 


LU 


0 


llll 


o< 


00 


□Q i— 4 


m<h 


O LL 


MQ^ 


LL« 


H 




1 — 


<0 




oO — | 


3LU 


0 


02 


LU LU 


LLO 


OO 


1 — 


LU 


H 


LUC 


00 


mmO 


2CQ 


3SO 


D_ M O 


CL — CD 


U-O 


H<< 


O 


00 


a 


H UJ 




X!50 


C 


o_j 


O*^ 


o>H 


c£UJa: 


00 


O 


< 


0 


00H 


Xm 


< 


LUO 


2 


Huj-r 


20 


<OCL 


ZQQ 


CL 


X 


CL 


200 


<x 


— LU LL 


>UJ 


CV) 


00OO 


X*“HLU 


HO 


hhLUlU 




CL 




*hZ 


< 


•^ooO 


CDQ 


COLU 


CC 2 


HO 2 


0O00 > 


CHH 


LU 




LU 


* — 1 


5: 


O •< 


P 


> 


M<cC 


O CL 2 


ocoarcG 


OHH 


CQ 


O 


00 


LU 


32 


hOIQ^ 


C r— S 


>M 


LL C 


MH C 


h-O-HLU 


<OQ 




2! 


< 


00 LU 


QCG 


UL •— 1 CL C 




CO CD 


LUX 


3:ox 


H Mh 


-JO 


O 


< 


X 


coo 


HX 


'—'LL LU • — 


• 00 00 




ULO 


o< 


00 • o'iu 


OCLCL 


H 




CL 


x< 


OH 


^OHO • 


OIZ 


2 — 


O 3 X 


2 2 LU 


mUDOO 


O LU 




LU 




CL X 


LUO 


ZlOHOHaJ< 


UJ 


OH 


<X cQ 


i-zo 


LUHLULU 




O 


O 


CL 


CLUJ 


O0 OO C <1 LL «H \ LL t-Q LU 


>00 


>20 


CQ O 


< UL LU 


CLHLUCD 


> 


2C 


2 


O 


i/jaOMM mwLL 


— 22 


♦“*X 


O »—( CL 


LUO 


H ^ cOqc*-hO 


H 


UJCZC 


2000 • 


O00GXXUJ w 


3 


OH 


2 LU 


OOl — 


oOUZ LU oO O O Q 


—4 


QCH 




<zz a 


L)oo<r<oQu 


U02O 


►O 


LU| — M 


LU 


O O O LO LU CL i — I — 


00 


UJUUJ 


cClu 


OO 


2UJ— »colu O 


IH 


3o0 


X<oO 


OC — t h-3Q 


X 


XlLO 


tu • — 1 1 — 


0 


ZZUJHXmxQ • 


H3I — 


■03 LU 


-« H 


QCGLU 


LU» 


LU 


oa 2 


O LU 1 — O 


2O00GGo:HcXOLUr-8 


CO 


LU2 C 


ULLU2 


LL 000; 




O 


OOO UJ 


ZUOI.U 


C Z LU CL! LU O < 


2H 


*—* 2 


X — 3 


O LU 


XhhLUcOO MM 


> 




OC 


UJ 2LUd: 


< 11 — 1 — X — 1>* 


►MOO lit 


<O0LL -J 


13m 


cohulx 


CL QOOOO 


JUHOLU 


CC-UCCcC 


X OOOOQ-OH 


LU 


CtcQH 


• O 


OOO 


:*>h <tt- h-i-^2: 


<z:z 


cL<rx 


LUCLlXO 


33Z 


ulLU<j 2U2DX 


2 U1MZ3 


:lloll 


U3I LU3 LULU 


Cd LU UJ 3 G 


xluqo 


cc 5 * — ' 


CL C.''U.U-'Zf-.00 


<JJlLlh M C 


l 2 LL 


O l-HL0O>QQ 


HG 


OOO 


axu 


hi 


OoO m3L11 


UJ2 


cQ ca H Cl 2 


:x-~*lu 


OH 1 O 


H 


oax 




OO LU 


OH O Q 33LOQ1M 


♦—4 


2 LU — 1 Cl — O! 



COGlUJ — IL-LIO*— 

<— JQ-CQUJuoCC 
<h^ H ^C^a: 
I — OOLLli-OH UJl — I — 



UJ LULU Ol— UUJ cQ 

cl 2 cl h 00 
OOLL r-* < UJ Q CL X < Q 
* It LM'JJhlh 



tuoo<Q cc 

CLUJ 0 < 0 <ClIiJXZIU 
^)air\ODza:u.a:<iua 
00 OOO LL LU 



1 — I — — I LU X — IcQOI — OUJ 
ZZJ(J *-h »-H 00 1— 

<OC zSZSDliD 
H2 LU *—* CO 5ZUJCL 



< 2*-< 


HO OOO OO HO X30 m (X 




cc 


LUjjL 


OOd >00oO 


Oo2H 0X2 




OoO 


O LULU UJHQoUJQX^ 00 


00 


o^coai 


crDtlO H 


ZQ.M< HC 


20 


02 — 


00 CL OCLCL dOO LU C)LU 


►— * 


LLO 


LL O UO CL LL O ZZ LU *-• 


OX 




CO 


• 30 


LL— <— 4 2 LU COOoO>U.QCLLL CvJO st200 


X 


LLO UL LU LLLU>2 


O LU LU*-i O CD 00 


MO 


OH ZXOd UwjMOH2n < 


< 


2 


— 12 






jXl 


00 •! — UL C 


^00«ZhOO2 C OOQTX 


1— * 


Hi 


m ( — iXOcD M IOX 


HCOhhLUXl CC 


ImmJ 


OHHLU HOOOOO oococl 




002-002 OOLU • 


OOOcLOUOoOLULLi- 


Oil- 0QU.O.UJOO3 00<<<o< a. 


LL 


m2 


mZmQmZCl(/)(3 


LL LL00CLLUC2 


CL X 1 — 1 0 


□ Q<C — 


hmCLQOQ 2 LU 0 LU 


►— * 


0 £*-*00 C£*-h 


UJ ♦— 1 LLI — • 




CL CL — iOOQ_ 


LU 


2 H 1 — 


OOO J LU ♦— * LU | — LU LUZOQ 




C 


HQ 


O IO O LU 


OO 


C oO UL *-h CL 2 1 CL 


3 oO C O —40 OC C2Zt/)ZlLZ<l LL 


> 


d x or oc a : llo • 


2 




C 


— 1 00 


2 O 


LLULULUL LLOC2CO<00 II 


O 


LL IU <1 UL LU 


• OLL OQ CL - 


2 II II 


II 




II 


II II 


II II 1 — 1 * — * * — • mCOOmoo uO t«l 


2 


OlL 


OQ-inooocL 0 


3 




II 




II H 


COUJ 


O 






X O *-u 


OH 


O 




2 


cC 00 


HrHCMcq >J- un no H lo h 




LUoO 


UJ 00 


OOUOH-HO 


OO 


or 


2 


0 


X<LUm 


0 11 11 11 11 11 11 11 11 c 


a 


3h 


3H 


UJHOOO 


OO 


UL 


O 


0 


CH IX 


0 X 


UJZJM 


— 1 H 


3*-»20 1 


2 


OC 


2 


2 


2 oOmC 


CL CL 


t-scz 


<C2 


02 LU O 


2 -* 


H 


C 


►— 4 


CD cQ LL *— « 


i — 1 1 — 1 


a:r>>3 


>3 


C 3 | O 


HH O 


00 


CD 




•—* mZ 




LUO 






> 2H 


*• 












oOO 






LU 


O 




in 




LOLOlA 




20 


O 


OO 


oO 0 X OC 


or 


in 


eg 




CsjrOsj- 




HH 


ocg 


covh 


LT\ sU H 3: UJ 


< 


H 


1 




1 1 1 ^ 


in m 


20 | 


1 V 


III LL 


0 


1 


vO 


O 


sj->4-vMn 


in 0 


Ot-H 


1 -« 


i — * »— # 


•— t rU t— 1 • UJ 


in 


O 


r—4 


ro 


(NJ rO sj- 




OC 


r— (rH 


c\Jco 


in 0 lQ oc 


H 












c 






0 


X 












0 






2 


LU 

2 









OOOOOOOOOOOOOO OOOOOOO^ OOO 000 OOOO OOOOOOOOOOOOOOOO 



122 



SCOR- 



OOOOOOOOOOOOOOOOOO OOOO OOOOOOOOOOOOOO 

00 O' O rH o J rT) UVO p- cO CX O ^ C\J ^ ^so^-co tro^OJinNt'-n O^OOCT^OrHCNj 

O'O'O'OOOOOOOOOQHHHHHt-trHHH r-HC\J C\J t\J k\l iNJ C\J C\J TnJ r\| (\J rH rH O 
O O r_ H rH H i— | rH r— 4 1 — I «-H H H •— 1 1 — 4 H H H ( H H i— I r—4 H <—4 >— 4 r— 4 r— 4 1—4 •— 4f— 4 H H pH 

OOOOOOOOOOOOOO OOOO OOOO OOOOOOOOOOOOOO 



ION 

MAX 

SS 

T 

ERE 

LIST 


2 

O 

t — • 




ED 

THE 

WILL 


22 


OOOO O oo X 


1 — 




►—I cO 


mUJ 


UJaiH c^LHOO^-* f— LU 


< 




G 00 < 2 


2 


1 — G 20 1 O 0 G X 


G 




CL 2 < O 




Z>Om >tO ►!- 


a 


G 2 


G q G m 


QG 


<21— oo gggoo 


CL 


00 G 


<h ! jjq «h— • 


GX 


2 < 2 oXO 2 • 


CL 


< X 


H* 0 1 — 2 O G 


001 - 


OOlG — • U l* — • I — 


G 


X 1 — 


GO DC 202 


00 


UJ 3 >Ci O G 00 


h- 


G 


CO G O Q < G Q 


GO 


GG mG£^ 002 X 2 cOm 


2 


•* 


G g gqc 


OH 


< 0 LL 0 cLUJu*-O 2 -J 


4—1 


x— a 


OcCLmCDOQ 


O 


GO G CD < G<Q 


— * 


O LU 


1 — QXg og 


CC G 


< -JC >2 001 — 12 G 




hXm 


g< clg 1 — 


GH 


o:ckolu ^3 *goo<x 


o<o 


XLUg 


OOGG Go 


< 


HGCL 2 (/) 20 <G !— 


CCGG 


2 X CL 


G < >LU 


G»— « 


OOO G 2 2 GGGOCXJ 


G 2 G 


l-CL 


stXUJ QhqI 


G DC 


G 3 C 0 <O'G 2 CL 02 2 


Nh'D 


1 - 3 


• 1 — X G — 1 i — Cl 


Cl 


CL XXOX— iO. <X»~i 


G 


<LLO 0 


O 1 — G 3 *— * O 


oo 


00 G 3 O h-<CH i — 


1 — O 


o 


rHOO OOOO 


H* CC 


G 2 G h- GtNJQCNJ 


U>G 


— * G • 


G 220 X 02 


Cl 


X 2 XG *COJ 0 X 


2 CD h— 


C 0 G 2 


1 — oo G 3 


LOG 


02 — 1 1 — GGQG 2 O 


00 


>r >^<2 


w< 


♦“•< 


— i<G*— <X • 2 O LUO 2 


UJQh 


• < 


LO| — G G — 1 < G 




1 X 23 mgZHUJZ< 


OOGG 


O O 00 G 


2 2 DC <2 LU X 


GO 


^02 oo G< 3 :< 


<2 


H| — G CO 


DO«Z •!— 


GO 


<<♦— • - Xooh- in 


X*— ' CO 


G HH 


hO 20 


»-<G 


D-OX 0 C GCJ 03 UJLH • 


cl<g 


Cl O i — 


1 — 00 oO <£ 2 Q 


G< 


o | ohsjh xcd «g 


4— 1 | ( — • 


— 3 ZL 


t — i — ZI<< 


O 


G< O 3 UJ XI— GO 


GO 


00 G G 


GGGOO—.G 


CL 


-GiXOOuOO 200 


LLOZ 


G "- *3 — 1 X G 00 00 *— • OG 


gg 



IU- CLt— I— M< #30 

UJQO(/)Um ujq CC 2 

2: vO uj x (XHZh 

ZLcC I 00dXDH<UM 
<GO'OOOOl— OL»-*GG 00 

icmno ♦-*!— gc\j 
os a: *:so>-oqi/) x 

3000 G — I — I h- 
<Zh GGG 2 G 00 X*-* 
XOCGooO 2001 - 3 : 
Li>OOMZ ZO H 
O g ro f-K/ >2 LL LU < QC 3 : Q 
i_ g | a : 

o g<o>: 50 x Oo Q< 
Z^oujDrg ogQ OoCU 
OGoo *oc uj < G 2 < 

ll cc 5: 3 m o cc LL UJ UJ LU < o cc 
H <3 C\JLL<OGCCX LU 
UZOO I G<J— C\J<X 

o »-«>fro omxc o h- 

X_J f\J — I | 3 X)<X — I LU O 

it < gqc -g>^gocqcggcd 2 
XZUJOZ .f~ hiz < 
UZCDOvJZoOUhOHZl- 



i — i LU t — • or O Q_ | — H- < Xo 
lu 3 oQLuiu<oixio<a:o h- 

>XO' Z'-UdCOlllOOItiJl H o LU 
G LU LUhhX <hdl“Z 2 0 

Z-J^ 3 — I | — OOG LU CC *— t LU <2 
O-jLL qclcoloqoo QO^ </) LU 
*-4 —LUQ-vt 3 o:< »zoOuj< 03 
(X 3 2 <XoO< LLcOX loo 

a: <OLLUO I GCL LU o cl lulu 

LUO G>- hh - 201 — goo 

00 LU O 2 -LU 3 02 — O I O OO Q 2 

2 *— ' G < CDG 2 2 LU LL 3 LU D 2 

— G*— XO LJLLOUJ 2*- 1 <Zl — LU 
Q ll — COO lu 3 QCODZ O O 22 

wQ_CL | H*XO | G< t 00 LU 2 
2 <G OOrH h-UJOO — ILUX 0022 CC < G 

o 3 oor-Lu QLf-r-ocoo occ gi— 

GoOOOG -G 2 LLg — LL GGt-.Q CD Ui • 
02 G • << oo • 2 <XI — O J 3 : 
Ooo • XH • LULU 01—0 LUCLq 
2 »— G 2 *OOf-LOX •>> LU UJ CD 51 — 1 
Oh -3 • 2 << 'ChZ -GX Og 

*— < O G — ' -00002 - 0 QO 3 — JCLh— 1-00 
*— G<o 00*-»00 GO< Q-GO 00 LL 



^<2 


1 < GGOOG <ooo 


GCZ > O C\1 1 — G G t— , CvJ <h-20LU 


3 < 


X 30 XCJG 3 CL 2 X 32 


GG CL | GG X | G 02 


00 


2 Q 


HO2 hOO 00 ZiOQOZ< 


GQ 2 <r-iG 


i 1 1 — 4 X 2 G <4. 


3 


GG 


2 


— 


GO GOgGG 


2 ^ 1 — O CL X 


< 


0003 


G 


m 


O G -dM 


G - LOOO 


O 


G<0 


3 


r —4 


OLU 3=200 


3 Oh- 




<OX 


arm 


| 


001— HU 


o-h< 




0 00 


G 1 


<r 


G< G | — | 


G 1 




X 


00 -T 


H 


OOlCDOOH 


OOrH 




1— OG 


G 




<G CC 


CD 




00 »-h 1 


3 




X HH 


3 




<IM 


00 




G G 


OO 




G^IL 



ooooo OOOO o OOOO OOOOOOOO OOOOOOOOOOOOOO 



123 



APPENDIX H 



ftftftftftftftftftftftftftftftft 

ftftftftftftftftftftftftftftftft 

ftftftftftftftftftftftftftftftft 2 



ft * ft 




ft ftft 




















•— § 


ft ft ft 




ft ftft 






















* * * 


oOUJ 


* * * 




















a 


ft ft ft 


UJCOO 


ft ft ft 




















UJ 


ft ft ft 


Z>OuO • 


* * * 




















CO 


ft ft ft 


-JCC-lh- 


ft ftft 




















3 


ft ft ft 


<CL<X 


ft ft ft 






















* * * 


> o 


ft ft & 




















oo 


ft * * 


•*_J t— « 


# ## 




















or 


ft ft ft 


ocO-jlu 


### 




















o 


ft ft ft 


<UUHl 


* 




















*- 


ft ft ft 


1— UJIS 


* * * 




















o 


ft ft ft 


00Q- 00 


* -> * 




















< 


ft ft ft 


3u023 


^ ** 




















u. 


ft ft ft 


<tO ' 


* 






















ft ft ft 


OOOCt* 


### 




















CD 


ft ft ft 


ZZOUJ 






















2 


ft ft ft 


<ImO> 


### 




















»— i 


ft ft ft 


:2 a: 






















-J 


ft * * 


2 CJL2 


* * * 


(NJ 


















< 


* * ft 


0-0 


* * # 


r-H 


















a 


■ft* ft ft 


^J>- LU _J 




w 


















oo 


ft ft ft 


i— • <OX ♦” * 


■R* "M* * 


Q 




















ft ft ft 


0021-00 


* 


LU 


















Q 


* * ft 


CLUJ Q_ 


* * * 


O 


















2 


ft ftft 


LU3 UJ 


* tt* 


H 


















< 


* * ft 


Of • 


^ -W- * 


t—4 




















ft ft ft 


LU'-UCtL 


* "tt* 


1— 


















oO 


ft ft ft 




* 


•» 


•# 
















O 


ft ft ft 


OU.C 




*— - 


oO 
















2 


ft ft ft 


CL Quo 


^ # # 


O 


1— 


















ft ft ft 


2 - h- 


* * * 


r-H 


2 
















H- 


* ft ft 


o>- 2 o 


** 




< 








CL 








1 — 


ft ft ft 


UhDJ 


***» 


UJ 


1 — 






X 


H- 




3 




UJ 


ft ft ft 


*— *hhCL 


* ft * 


—1 


00 


O 




o 


O 




< 




UO 


* ft ft 


OoO 1— 


* * * 


h- 


2 


>J- 






—J 




CC 






* ft ft 


H-2<Q- 


* 


»— i 


O 


• 




• 


CL 




Q 




2 


* * ft 


luclx 


* * * 


h- 


O 


o 




o 












ft ft ft 


QOCQQ 


-R- ** 


•» 








•* 


UJ 




UJ 




< 


* ft ft 


LU ♦— i O 


^ R- # 


N 


o 


II 




00 


2 




2 




O 


ft ft ft 


2JJJ 


* * # 


m 


2 






o 


tH 




►— « 






ft ft ft 


cd<<< 






»— t 


1 — 




• 4 * 


H* 




H* 




uO 


ft ft ft 


»-»oroo 


# ft# 






2 




o 


Z> 




Z) 




3 


ft * ft 


oO|— 


X <- * 




o 


C 




• 


o 




o 




O 


* * ft 


UJO>Q 


*«• * * 




—1 


1— 




CNJ 


QC 




a: 






ft * ft 


OUJh-2 


## ^ 




-J 


00 




X 


OQ 




CQ 


o 


cC 


* * ft 


CLhhC 


* * * 


w 


o 


2 


no 


CC CO 


z> 


LU 


r> 


UJ 


< 


ft ft ft 


OOOOo 


# ^ * 


> 


X 


O 


O 


< ft 


oO 


—i 


00 


-J 


> 


ft * ft 


<t Do: 


# * * 


•» 




o 


>t 


^ft 0> 




I— 




H* 




ftft ft 


I2LL_JLU 


* ^ * 




LU 




O 


2oo in 


CC 


*— « 


CL 


»— i 


UJ 


* * * 


OUJI— 


# * * 




X 


2 


• 


O'J' H 


a 


i — * 


O 


H* 


X 


* * * 


2 >2 


*** 


-Jw 


1— 


< 


CJ 


r— ( nJ" 


X 


00 


U- 




1 — 


* * ft 


<t oO »-* 


^ -w- * 


LUX 




2 




— • r— i ro 




— < 




-^ao 




* * ft 


or h-oot 


* ^ * 


CD 


UJ 


or 


II 


<o • ► 


UJ 


OO 


UJ 


in< 


UJ 


ft ft ft 


0 02 0- 


^ * * 


<Z 


M 


< 




X rOH 


_J 


fO ■ — 4 


-J 


nOO 


M 


ft ft ft 


OM 


* * * 


-*o 


v— * 




< 


CL II ft II Q O UJ 


1 — 


»*w 


H- 




►— < 


* * * 


aZ uj 


R* ^ # 


r-H 


_J 




X 


— 1 ••— • • *Q 


t— « 


in 


»— i 


in 


-J 


ft ft ft 


Q-^|— Q 


#*# 


00 00 


< 


2 


CL 


<QN OO^ 


1— 


'^i — 


H* 


^i— 


c 


ft ft ft 


X*-* 




ft 2 


>— < 


a 


_J 


LU O II II ^ 




< 




< 


1— « 


* * * 


oozo> 


# 


— J LU 


H- 


> 


< 


<co ii r\|~.^|_ 


Q 


Q2 


Q 


02 


1 — 


ft ft ft- 


♦^i O ♦— «o 


* * * 


<12 


►— » 






l-Z) mmZ 


< 


<icc 


<i 


cor. 


1 — < 


ft ft ft 


Xcxluo: 


# * * 


LU^-h 


2 






<iumo^^o 


UJ 


LUQ 


LU 


LUO 


2 


ft ft ft 


J-LLXQ- 


### 


o:q 








Q>XQX>0 


or 


aiLL 


CL 


oru. 


♦— i 


* * * 




* * * 






















ft ft ft 




* * * 










O 




o 




in 




ft ft ft 


ft ft ft ft ft ft ft ft ft 


ft * * * 










04 




m 




CO 





ft * ft * ft ft ft ft * ft * ft ft * ft ft 
ft * ft ft ft ft ft ft ^ ft ft ft ft ft ft ft 
oooooocjcjooooooooocjo 



OOOOOOOO 



OOO 



OOO 



uo 



2 

o 

3 

CO 

CC 

UJ 

> 

2: 

a 

o 

CO 

o 3 

Z C 

*-« Q h- 



3 

CL 

K— 

3 

O 

LU 

o 

< 



o 

2: 

UJ 

> 



o 

UJ 

UJ 

CL 

to a: 

a 

03 

2:0 

*-*< 

J3CLL 



3 


QC 


»— 1 


O 




< 


CO 














3 


O 


0 


>3 




3 


UJ 3 














UJ 


O 


►— 1 


CC 




w 


X3 














to 


UJ 


0 


X3 




•* 


< 
















a: 




o> 




3 


►a: 














UJ 




0 


oc2 




3 


>-3 




O 










X 


O 


1 — 


LLO 




3 


03 




2 






to 




3 


3 




O 




•* 


23 




3 






CD 








e> 


3 




m 


3< 




3 






O 




LL 


cc 


0 


QCO 




XI 


30 




< 






CO 




O 


O 


3 


LU vO 




• 


O 




O 












3 


< 






•* 


LU> 




to 






3 




2: 


a: 


2 : 


2 : 




3 


OC 3 














O 


CL 


< 


ox 




LA 


33 




3 






DC 




1— 1 






oca 




1 — • 


O 




< 




to 


< 




3 


2 : 


0 


i — 1 






►O 




2 




3 


H 




CL 


LJ 


z 


0 




m •*"* 


X. 3 




a 




O 


to 




1— 1 


i — 1 


1— 1 


30 




QCtO 


33 


to 


3 




3 


3 




CC 


H 


a: 


3 






CL> 


3 


CO 




< 


9* 




0 


< 


3 


Q 




CQ2 


3 


3 


2 




> 


C\ l 




to 


O 


a 


•30 




30 


00 


CO 


3 






tn 0 




LU 


i — 1 




COO 




• 3 


LU2 


3 


X 




CZ 


• 3 




O 


LL 


2 : 


3rO 




► 3 


CL<t 


•» 


KH 




< 


3 to 






i — 1 


0 


O' 




O' CO 


to 


O 


a 






# CL 




< 


3 


3 


CC 




roa. 


to 








CO 


# 3 






CL 


3 


Oco 




3LU 


>-3 


2rO 


3 




0 


— • ► 






X 


C 


30 






03 


► • 


x 






r- 0 




• 


< 


O 


OcO 




m •» 


lZ> 


XlO 


1 — 




a 


0 




3 




1— 1 


< 




00 


33 


► 3 






2 


O rO 3 — 




2 


3 


UL 


33 




304 


23 


CC 3 


CO 




< 


0 no >r- 




LU 


< 


3 


X 




a: 3 


3 


< ► 


3 






0 co < • 




X 


2: 


3 


UJH 




33 


3 


CQ<0 


a 




•> 


0 CO 20 




1— 1 


O 


CL 


3 




m ► 


3G0 


3 • 


»— < 


rO 


2 


• rO »>3 




cC 


►— 1 


V 


CO 






xo 


►O 


> 




0 


O • 5X1 3 




LU 


CO 


< 


02 




4- a: 


y-cc 


33 


0 




3 


* 0 a ► ► 




CL 






to^ 




(XI < 


Cl 


►3 


a : 


— * 


3 


* # 3CCX 




X 


UJ 


LU 


CC 




33 


CC 


3 *» 


3 


O 


CO 


0 ■»■ cocro 




UJ 


X 


X 


33 




►to 


03 


3 3 






3 


2 3 30 Qw 






3 


3 


IQ 




- 3 


30 


M 1 — < 


CC 


<xj 


3 


3 M oox- 




LU 






3 




3 - 




X ► 


0 






> ► ► 




X 


to 


to 


3 




3 ► 


too: 


Q-rO 


1 — 


t—i 


•* 


< QC<t-^3>t 




3 


3 


3 


C03 




3>t 


33 


• 


0 


< 


a: 


3 CcO— •• • — 






2 : 


X 


3 cl 


— * 


IO 


3cO 


— O 


< 


O 


3 


XO>t33 (NJ 




LL 


UJ 


LU 


Zh 


rO 


CL 3 


3X 


O 3 


LL 


* 


CO 


<X « j 2 X ^-33 O 




O 


CO 


CO 


33 


nj 


• 3 *— 


<3 


03 




3 


r 


CO X<OCO^HHLU 3 






UJ 


LU 


LOO 


* 


► - X 


>2 


It ro 


2 


2 


O 


O 3O>rj0t/)X ► tO 




to 


CC 


CC 


3 


* 


CO- 0 




O ► 


U 


t—i 


2 


\ ^ K & Jcaam^a 




LU 


CL 


CL 


<X3 


Osi 


32 3 


C ^ 


20 


1 — * 


< 




3 <3(NC03 3 M 3 




O 


UJ 


LU 


D-< 




*— 1 » ►O "'*• 


33 


3 • 


3 


0 


3 


^ xa03co-» 


0 


< 


CC 


OC 


LU3 


X 


O- 30- O 


<x 


►O 


< 


* 


> 


H3dJJ^ ^O- OO 


0 


3 






OCe-f 


(XI • 


vt3Minom 


QO 


O 3 


or 


3 


< 


3 3 J CO tO 3 •-'OO O 3 3 


*3 


to 


H 


(XI 


O 


300 


•**0 CO ►— •* 


*— • 


X 3 


cO 


3 




m<q.q.cl Oo 


•> 


• • 






3h3 


• *o 


vO^CL -O'— 0 


3 3 


•-w 


3 


< 




llX^-33 ll vO^OU 


0-0 


to 00 


2 


2 


3QO 


003 


— » UJ ' 


XX 


in 


3 


O 


3 


3 II II ^ II ^ — 1 3 


4)w 


3 3Z 


►— 1 


M 


< C 


1 


3- 3 


3 


~3 


< 


CO 


1 — 


O-tt 3 (XI II tO h-<< 




OO 


< 


< 


UOh 


3c\]3 


3C -3<3 


3 


< 


O 


■w 


O 


2300 3 coca 3 < II H 303 


Hj 


0 


O 


tOH CO 


ZZJ 


h^COH^h 


03 


ax 




II 


3 


3 11 JJ<cQDh2:<-^ + 


1 — i — 


OT3 








MH «C 1 


3 3 cC CO 3 Q1 3 


CO 


< CC 


3 


3 


X 


<tO</)l — OcO»— 4 OC*~‘»— I»— < 


t—i 


<0 








<<o h 


3QC 


30 


X 


< 


u 


<IQ.a^lO^CCO w '- II Orz 


>LL 








OOoOm^LL ►2llz 


OCCL 


CC 3 


h 


O 


0 


3:33 3 O JC3 QC 3X >3 03: 












1 

O O 




OO 








O 


0 












vf in 




03 








CO 





UOOOOOOOOOO 



0000 000 000 



125 



</> 

< 



o 

•* 

cr» 

•* 

o 

9» 

o 

•* 

o 

•* 

o 



o 









•* 




• 




in 




00 




• 




UJ 




0 




O 




•* 




_J 




in 




< 




• 




> 




0 


• 






•* 


00 


> 




0 


IU 


m 


— > 


UJ 


I- 


•» 


r*- 


_j 


< 


00 


• 


h- 


2 


rO 


m 




t— 1 


h- 


•h 


1— 


Q 


•» 


LL 


•* 


r£ 


. 




_J 


a 


oo 


X 


LU 


0 


LU 


0 


lQ 


0 


O 


■— * r— 1 


< 




_J 


*— 1 •* 


_J 


a 


< 


wf^- 


•* 


►— * 


> 


>- • 


— rH 


a: 




-cn 


O LU * 


0 


X 


— % rH 


O 


• 


. 


H^U- 




•» 


•» 


•w •» 




0 


in 


xo 


x>-^ <t . 


CNJ 


rH 


r— < 


0 - ox 


1- 


*"*1— 






**210 •* 


O *“*— -o •‘PO 


- 0—1. 


*.f\J . 


> 4 - • m- ^ 


HlA^O'HrHO 


*H Q, rH O 


*» 


•* 


II .. 


f— 5T 



O > 0 '-'^ 'O— 'UJsO^O^O^-C 



— — — r>*- -J^ ^ 

h- h- o z: j~CL i — q 

<LULU<fOUJ<*— UJ< LU< 

s: j- i- s: rH h- ^ h- i- a : -j i— 2 : -j a. 

a I mmix; hc/zmqcjmo: joq 

o q: oc o o cc o o cc q <r a. a <1 h* 2 

LL 3: 3: UL Q is LL O 3 : UL O 12 LL O 00 LU 



O 


O 


OO 


O 


O 


O 


rH 


rvjrn 


<r 


m 


rH 


rH 


rHrH 


rH 


•H 



126 



LIST OF REFERENCES 



1. Busch, N. E. , 1973: The surface boundary layer. 

Bndy-Layer Meteor . , ^4, 213-240 . 

2. Cavanaugh, M c P., 1974; Examination of shipboard 

measurements of the vertical profiles of mean tem- 
perature, humidity and wind speed. M.S. Thesis, 

U. S. Naval Postgraduate School, Monterey, California. 

3. Davidson, K. L., 1970; An investigation of the in- 
fluence of water waves on the adjacent airflow. 

Ph.D. Dissertation, University of Michigan, 259 pp. 

4. Davidson, K. L., 1974: Observational results on the 

influence of stability and wind-wave coupling on 
momentum transfer and turbulent fluctuations over 
ocean waves. Bndy-Layer Meteor . , (in press), _5, 123- 

145. — — 

5 . Davidson, K. L. and Frank, A. J., 1973: Wave-related 

fluctuations in the airflow above natural waves. 

J. Phys. Oceanogr ., 3_ ( 1 ) , 102-119. 

6. Garratt, J. R. , 1972: Studies of turbulence in the 

surface layer over water (Lough Neagh). Part II. 
Production and dissipation of velocity and tempera- 
ture fluctuations. Quart. J., Roy. Meteor. Soc ., 98 , 

6 4 2-657. " — - 

7. Haltiner, G. J. and Martin, F. L., 1957: Dynamical 

and Physical Meteorology , McGraw-Hill Book Company, 
(New York), 470 pp . 

8. Kaimal, J. C., Wyngaard, J. C., Izumi, Y., and Cote, 

O. R. , 1972: Spectral characteristics of surface- 

layer turbulence. Quart. J., Roy. Meteor. Soc ., 98 , 
563-589 . 

9. Kraus, E. B., 1972: Atmosphere-Ocean Interaction , 

Clarendon Press, (Oxford), 275 pp. 

10. Lumley, J. L., and Panofsky, H. A., 1964: The 

Structure of Atmospheric Turbulence , I n t er s c i enc e , 

(New York), 239 pp . 



127 



11 . McKendrick, J. M . , 1972: An investigation of digital 

spectral analysis programs and computer methods 
utilized at the Naval Postgraduate School in the 
analysis of high frequency random signals. M.S. 
Thesis, U. S. Naval Postgraduate School, Monterey, 
California, 147 pp . 

12 . Paulson, C, A., 1970: The mathematical representa- 

tion of wind and temperature profiles in the unstable 
surface layer. J. Appl. Meteorol . , 9_, 857-861. 

13 . Pond, S., Smith, S. D. , Hamblin, P. F. and Burling, 

R. W. , 1966: Spectra of velocity and temperature 

fluctuations in the atmospheric boundary layer over 
the sea. J. Atmospheric Sci ., 23_(4), 376-386. 

14. Stegen, G. R. , Gibson, C. H., and Friehe, C. A., 

1973: Measurements of momentum and sensible heat 

fluxes over the open ocean. J. Phys . Oceanog ., _3(1), 
86-92. — 

15. Taylor, G. I., 1938: The spectrum of turbulence. 

Proc o Roy . Soc . , A164 , 476 . 

16. Welsh, P. T., 1974: An investigation of ship related 

motions on turbulent measurements. M.S. Thesis, 

U. S. Naval Postgraduate School, Monterey, California. 



128 



INITIAL DISTRIBUTION LIST 



No • 



1. Defense Documentation Center 2 

Cameron Station 

Alexandria, Virginia 22314 

2. Library (Code 0212) 2 

Naval Postgraduate School 

Monterey, California 93940 

3, Naval Oceanographic Office 1 

Library (Code 3330) 

Washington, D . C. 20373 

4, Commander, Naval Weather Service Command 1 

Naval Weather Service Headquarters 
Washington Navy Yard 

Washington, D . C. 20390 

5, Professor Kenneth L 0 Davidson, Code SIDs 20 

Department of Meteorology 

Naval Postgraduate School 
Monterey, California 93940 

6. Professor Thomas M. Houlihan, Code 59Hm 3 

Department of Mechanical Engineering 

Naval Postgraduate School 
Monterey, California 93940 



7. Mr. P. Vial, Code 048 
Naval Ordnance Laboratory 
White Oak 

Silver Spring, Maryland 20910 

8. Mr. E. Boudreaux, Code 048 
Naval Ordnance Laboratory 
White Oak 

Silver Spring, Maryland 20910 

9. Dr. P. Livingston 

Applied Optics Branch, Bldg. 30 
Naval Research Laboratory 
Washington, D. C. 20390 

10. LT William E. Johnston 
c/o Mr. Lou Hanna 
33 East Bond Street 
Corry, Pennsylvania 16407 



1 



1 



1 



3 



Copies 



129 



1 



11. Mr. Steve Rinard 
Department of Meteorology 
Naval Postgraduate School 
Monterey, California 93940 

12. Mr. Robert Smith 

Department of Research, Code 023 
Naval Postgraduate School 
Monterey, California 93940 

13. Ms. Sharon Raney 
Computer Center, Code 0211 
Naval Postgraduate School 
Monterey, California 93940 

14. Mr. Robert Limes 
Code 7 2Ec 

Naval Postgraduate School 
Monterey, California 93940 

15. Professor Dale F. Leipper 
Chairman, Department of Oceanography 
Naval Postgraduate School 
Monterey, California 93940 

16. LT Michael P. Cavanaugh 
c/o Mr. James A. Cavanaugh 
3607 Husted Drive 

Chevy Chase, Maryland 20015 

17. LT Patrick T. Welsh 

Student, Naval Destroyer School 
Naval Station 

Newport, Rhode Island 08240 

18. R/V ACANIA 

Department of Oceanography 
Naval Postgraduate School 
Monterey, California 93940 



1 



1 



1 



1 



1 



1 



1 



130 



SECURITY CLASSIFICATION OF THIS PAGE (HT,™ Entered) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. REPORT NUMBER 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG NUMBER 


4. TITLE (end Subtitle) 

Estimating Boundary Layer Fluxes fr 
Measurements of Dissipations of Tur 
Energy and Temperature Variance 


om Shipboard 
bulent Kinetic 


5. TYPE OF REPORT * PERIOD COVERED 

M.S. Thesis 
March 1974 




6. PERFORMING ORG. REPORT NUMBER 


7. author^; 

William Edward Johnston 


8 . CONTRACT OR GRANT NUMBERS) 


9. PERFORMING ORGANIZATION NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 93940 


i 10. PROGRAM ELEMENT, PROJECT, TASK 
AREA & WORK UNIT NUMBERS 


II. CONTROLLING OFFICE NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 93940 


12. report DATE 

March 1974 ' 


13. NUMBER OF PAGES 

132 


t4. MONITORING AGENCY NAME & ADDRESS^// different from Controlling Office) 

Naval Postgraduate School 
Monterey, California 93940 


15. SECURITY CLASS, (of thie report) 

Unclassified 


15 a, DECLASSIFICATION/DOWNGRADING 
SCHEDULE 


16. DISTRIBUTION STATEMENT (of thie Report) 

Approved for public release; distribution unlimited. 


17. DISTRIBUTION STATEMENT (of the abatract entered In Block 20, It different from Report) 


18. SUPPLEMENTARY NOTES 


19. KEY WORDS (Continue on reveree aide If neceeeery and Identify by block number) 

Turbulence Friction Velocity 

Marine Boundary Layer 

Dissipation of Turbulent Kinetic Energy 
Momentum Flux 


20. ABSTRACT (Continue on reveree aide It neceeeery and Identify by block number) 

i Velocity and temperature fluctuation measurements were made over the open 

ocean from instruments mounted on the R/V ACANIA. These data were examined 
to determine the validity of present formulations and prediction techniques. 

Values of momentum flux, u* 2 , were inferred from the rate of dissipation 
of turbulent kinetic energy, e. Dissipation values were obtained from spectral 
estimates and inner scale estimates. Values of u* were examined for repre- 
sentativeness on the basis of the constant-flux assumption and by comparisons 



DD i jan^73 1 473 EDITION OF l NOV 65 IS OBSOLETE 

(Page 1) S/N 0102-014- 6601 | SECURITY CLASSIFICATION OF THIS PAGE (»T<.n D»(» Bnl.r.d) 

131 




SVCUHITY CLASSIFICATION OF THIS PAG EfWTien Dmtm Entmrmd) 



with other studies. The vertical variation of the dissipation rate was 
examined for possible effects of stability and wind- wave coupling. 

The momentum flux, computed from spectra, supported the constant-flux 
assumption for neutral conditions. For periods of instability, the 
stability corrections applied to the vertical variation of £ resulted in 
the proper adjustments toward the predicted slope. The reductions in 
momentum transfer during periods of stable stratifications were consistent 
with wind-wave coupling effects described by Davidson. The shape of a 
spectrum of temperature fluctuations was in agreement with predictions. 



DD Form 1473 
, 1 Jan 73 

S/N 0102-014-GG01 



132 



(BACK) 



SECURITY CLASSIFICATION OF this P AGEfWTi.n O.l. Enffd) 



r 



Thes i s 

J675 

c.l 



23 SEP 74 



2 110 7 



Thesis X^' r>f ^PQ 

J675 Johnston 

C,:L Estimating boundary 

layer fluxes from ship- 
board measurements of 
dissipations of turbu- 
lent kinetic energy and 
temperature variance. 

23 SEP 74 '21107 



11*0586 

Johnston 

Estimating boundary 
layer fluxes from ship- 
board measurements of 
dissipations of turbu- 
lent kinetic energy and 
temperature variance. 



